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ABSTRACT 


Fluctuations  in  the  amplitudes  and  phases  of  20-sec  Rayleigh  waves  were 
measured  for  ten  earthquakes  recorded  at  NORSAR  and  for  two  earthquakes 
recorded  at  an  1100-km  long  linear  array  in  the  Southwestern  United  States. 
These  measurements  were  compared  with  the  prediction  of  Chernov  (1962),  who 
analyzed  the  scattering  of  elastic  waves  by  random  inhomogeneities  along  the 
source-to-receiver  path.  Use  of  this  theory  enabled  the  amplitude  and  phase 
fluctuations  to  be  related  to  a  statistical  description  of  the  random  scatter¬ 
ing  medium.  As  predicted,  it  was  found  that  the  fluctuations  were  correlated 
over  a  longer  distance  in  the  direction  parallel  to  the  propagation  of  the 
wavefront  than  in  the  direction  perpendicular  to  it.  This  implies  by  recipro¬ 
city  that  relative  Mg  within  a  test  site  would  be  better  determined  at  a 
station  in  line  with  the  vector  between  two  events  than  by  a  station  perpendi¬ 
cular  to  the  vector.  Certain  measurements  were  incompatible  with  Chernov's 
theory,  however;  these  discrepancies  may  be  attributable  to  multipath  arrivals 
and/or  large-amplitude  scattering  for  which  the  Born  approximation  is  not 
valid. 

It  was  found  that  narrow-band  spectral  measurements  of  Rayleigh  wave 
amplitudes  exhibited  stronger  fluctuations  across  both  arrays  than  did  visual 
measurements  of  Mg.  The  size  of  the  fluctuations  was  diminished  by  using 
broad-band  spectra.  Fluctuations  in  Mg  are  approximately  1/3  those  in 
across  equivalent  array  dimensions.  By  reciprocity  this  implies  that  relative 
yields  within  a  test  site  may  be  determined  by  using  Mg  with  a  standard  devia¬ 
tion  1/3  of  that  obtained  from  m^  if  an  equivalent  number  of  station  measure¬ 
ments  are  available. 

Random  scattering  of  Rayleigh  waves  has  a  larger  effect  upon  narrow-band 
measurements  of  attenuation  than  upon  measurements  of  phase  velocity.  As  a 
result  of  such  scattering,  surface-wave  spectra  may  display  significant  fluctu 
ations  at  various  frequencies,  causing  source  mechanisms  inferred  from  spectra 
observations  at  only  a  few  stations  to  be  unreliable. 
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INTRODUCTION 


Researchers  have  long  recognized  the  effects  of  the  earth's  lateral  inho¬ 
mogeneity  on  basic  seismological  measurements.  In  fact,  systematic  variations 
of  wave  amplitudes,  travel  times,  and  spectral  content  have  been  successfully 
utilized  to  infer  much  of  the  structural  variations  of  the  earth's  crust  and 
mantle.  However,  in  most  cases  the  structure  becomes  discernible  only  when  a 
suitable  amount  of  data  is  averaged,  implying  that  a  significant  degree  of 
randomness  appears  in  that  data  and  that  single  observations  are  not  reliable 
for  making  geophysical  inferences.  For  example,  the  range  of  P-wave  amplitudes 
over  LASA  is  frequently  as  large  as  10  to  1  (Chang  and  von  Seggern,  1977), 
with  the  amplitude  pattern  changing  among  events.  The  random  component  in 
data  sets,  such  as  travel  times  or  amplitudes,  can  usually  be  described  with 
common  statistical  formulations  and  a  statistical  description  of  the  media 
properties  can  often  be  derived  from  the  statistical  content  of  observed  seis¬ 
mic  data.  This  observed  random  scattering  is  associated  with  features  of  the 
medium  having  smaller  scale  than  those  causing  broad  regional  effects,  such 
as  North's  (1978)  magnitude  bias  or  Herrin  and  Taggart's  (1968)  travel  time 
residuals.  However,  because  of  its  severity,  random  scattering  often  obscures 
broad  regional  effects  and  inhibits  determining  underlying  structure  and  seis¬ 
mic  source  parameters. 

Aki  (1973)  and  Capon  (1974)  performed  scattering  analysis  on  the  travel 


Chang,  A.  C.  and  D.  H.  von  Seggern  (1977).  A  study  of  amplitude  variations 
and  mb  bias  at  LASA  subarrays,  SDAC-TR-77-11 ,  Teledyne  Geotech,  Alexan¬ 
dria,  Virginia. 

North,  R.  G.  (1978).  Station  bias  for  ISM-reported  magnitudes,  Geophys.  J., 
in  press. 

Herrin,  E.  and  J.  Taggart  (1968).  Regional  variations  in  P  travel  times. 
Bull.  Seism.  Soc.  Am. ,  58,  1325-1337. 

Aki,  K.  (1973).  Scattering  of  P  waves  under  the  Montana  LASA,  J,  Geophys. 
Res.,  78,  1334-1347. 

Capon,  J.  (1974).  Characterization  of  crust  and  upper  mantle  structure  under 
LASA  as  a  random  medium,  Bull.  Seism.  Soc.  Am.,  64,  235-266. 


times  and  amplitudes  of  LASA  P-wave  recordings  using  Chernov's  (1962)  theory. 
This  report  focusses  on  whether  or  not  a  similar  phenomenon  exists  for  Rayleigh 
waves  and  how  it  can  be  described  theoretically;  it  will  not  deal  with  the 
phenomenon  of  distinct  multipathing  of  surface  waves  that  Capon  (1970)  studied 
because  multipathing  is  a  deterministic,  not  statistical,  effect  arising  from 
known  lateral  variations  in  the  gross  structure  of  the  crust  and  upper  mantle 
such  as  the  oceanic/ continental  boundary.  However,  such  multipathing  should 
grade  into  random  effects  on  surface-wave  amplitudes  and  phases  where  the 
structure  is  uniform  except  for  small-scale  heterogeneities  that  can  be  des¬ 
cribed  statistically  using  Chernov's  method.  Demonstrable  random  scattering 
effects  on  Rayleigh  waves  will  have  important  implications  for  the  reliability 
of  surface-wave  magnitude  estimates,  for  the  accuracy  of  surface  wave  phase- 
velocity  and  attenuation  measurements,  and  for  the  character  of  surface  wave 
spectra. 

The  first  portion  of  this  study  examines  evidence  that  suggests  random 
scattering  of  Rayleigh  waves,  then  modifies  Chernov's  three-dimensional  theory 
to  predict  the  effects  for  a  two-dimensional  propagation  problem  and  presents 
new  results  that  can  be  compared  with  the  modified  theory  and  other  possible 
models.  Finally,  the  probable  error  of  various  Rayleigh-wave  measurements 
taken  from  long-period  seismograms  is  investigated.  Throughout,  waves  with 
periods  near  20  sec  are  stressed  because  cf  the  practical  importance  of  using 
surface-wave  magnitude  (Mg)  in  estimating  yields  of  underground  nuclear 
explosions . 


Chernov,  L.  A.  (1962).  Wave  propagation  in  a  random  medium.  New  York,  Dover 
Publications,  Inc. 

Capon,  J.  (1970).  Analysis  of  Rayleigh-wave  multipath  propagation  at  LASA, 
Bull.  Seism.  Soc.  Am.,  60,  1701-1731. 


OBSERVATIONS  OF  SMALL-SCALE  SCATTERING  OF  RAYLEIGH  WAVES 


Few  opportunities  exist  to  observe  recorded  surface  wave  signal  variations 
over  small  spatial  dimensions  on  the  order  of  a  few  wave-lengths  or  less  with 
a  regularly  spaced  array.  Aside  from  the  large  seismic  arrays,  LASA  and  NOR- 
SAR,  there  are  almost  no  available  data  sources.  However,  during  early  1962 
one  well-placed  set  of  stations  existed  in  the  Southwestern  United  States. 

These  sites,  part  of  the  VELA  Long  Range  Seismic  Measurements  program  (Figure 
1),  were  situated  to  record  seismic  signals  from  the  GNOME  underground  nuclear 
explosion  and,  all  together,  composed  a  1100  km  linear  array.  During  their 
short  recording  interval,  surface  waves  from  several  earthquakes  at  teleseismic 
distance  were  well  recorded  at  most  of  the  11  sites.  For  this  study,  signals 
arriving  perpendicular  and  parallel  to  this  linear  array  were  needed  to  verify 
certain  theoretical  models;  two  such  events  were  available.  Their  epicenter 
information  and  great  circle  paths  are  indicated  on  Figure  1;  Figure  2  and  3 
illustrate  the  vertical-component  recordings  of  the  Rayleigh  waves  for  the  two 
events.  Note  that  the  plots  preserve  the  true  relative  amplitudes.  Signifi¬ 
cantly,  the  range  of  amplitude  among  the  seismograms  at  a  period  of  20  sec  is 
over  a  factor  of  three  and  it  is  apparently  random.  Attenuation  along  the  path 
or  changes  in  crustal  responses  at  the  individual  sites  cannot  alone  explain 
this  range  and  a  detailed  analysis  of  these  particular  recordings  will  be  pre¬ 
sented  in  this  report. 

Observations  of  random  fluctuation  in  Rayleigh  waves  are  possible  at  the 

NORSAR  array  of  long-period  instruments,  (shown  in  Figure  4),  but  only  over  a 

span  of  100  km.  One  example  is  the  recordings  made  of  the  Turkish  earthquake 

shown  in  Figure  5.  Comparing  the  seismograms  with  those  of  Figures  2  and  3 

reveals  less  variation  in  signal  amplitudes  and  shapes,  an  observation  to  be 

expected  because  the  aperture  is  reduced  by  an  order  of  magnitude.  Yet,  among 

the  NORSAR  recordings  there  are  significant  differences  that  cannot  be  explained 

by  varying  crustal  response  or  slight  calibration  errors.  For  such  typical 

events,  the  standard  deviation  of  measured  M  over  22  NORSAR  sensors  will  be 

s 

somewhat  less  than  0.1.  A  thorough  analysis  of  these  types  of  recordings  will 
be  presented  later  when  we  compare  theoretical  and  empirical  results  for  the 
random  scattering  of  Rayleigh  waves. 
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Figure  3.  LRSM  long-period  recordings  for  the  Kermadec  earthquake. 


NORSAR  LP  sensor  effects 
.0508  hz  (1  smoothing) 
10  events 


10  km 


Figure  4.  Layout  of  the  NORSAR  long-period  array.  Numbers  indicate  sensor 

effects,  the  difference  between  the  individual  spectral  amplitudes 
and  the  mean  spectral  amplitude,  averaged  over  10  events. 
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Rayleigh  waves  recorded  at  NORSAR  from  an  earthquake  in  Turkey  (06 


Using  wavenumher  analysis  Mack  (1972)  has  already  investigated  this  ap¬ 
parent  scattering  phenomenon  at  arrays,  suggesting  that  actual  Rayleigh  wave 
signals  can  be  represented  in  two-dimensional  wavenumber  space  as  a  distribu¬ 
tion  rather  than  as  a  single  discrete  point.  One  of  these  possible  representa¬ 
tions  is  the  Gaussian  distribution: 

F(kx)  =  exp  (-a2kx2);  F(ky)  =  exp  [~S2(ky  -kQ2]  (1) 


where  the  kx  and  k  axes  are  parallel  and  normal  to  the  wavefronts,  respectively, 

and  k  is  the  wavenumber  where  the  power  of  the  signal  peaks  for  a  given  fre- 
o 

quency.  From  this  model.  Mack  derived  the  form  of  the  coherence  function  for 

spectral  amplitudes  from  two  sensors  at  the  frequency  corresponding  to  k  as 

2  2  2  2  ° 

Y2(r")  =  [exp  -  V  rx  .  exp  -  W  ry]  (2) 

a2 

wiiere  r  is  the  vector  from  one  sensor  to  another  within  the  array  of  instru¬ 


ments  and 


r  =  (r  ,r  ) 
x*  y 

rx  =  X1  _X2 
ry  =  yl  "y2  ‘ 


(3) 


To  simplify,  henceforth  r  =  jr|.  At  the  frequency  of  .047  Hz,  Mack  computed 
coherence  between  sensor  pairs  for  the  signals  shown  in  Figure  5;  the  results 
are  in  Figure  6.  Results  for  sensor  pairs  nearly  normal  or  parallel  to  the 
wavefront  are  shown  as  a  function  of  sensor  separation  r.  We  have  superimposed 
the  relation  (2)  for  various  of  a  (or  3).  Assuming  rx  =  0  for  sensor  pairs 
aligned  normal  to  the  wavefront,  we  estimate  3  ^  1200;  similarly  a  'v  500  for 
sensor  pairs  aligned  parallel  to  it.  Clearly,  coherence  along  the  raypaths 
falls  off  much  less  rapidly,  a  fundamental  feature  of  the  random  forward  scat¬ 
tering  of  waves  detailed  in  Chernov's  (L962)  theory.  However,  alternative 
wavenumber  representations  or  signal  models  will  predict  similar  coherence 
patterns;  for  instance,  two  discrete  multipath  arrivals  of  nearly  equal  ampli¬ 
tude  and  equal  phase  delay  at  nearby  points  k^  and  k£  in  wavenumber  space 
would  result  in  a  theoretical  coherence  function  (Mack,  1972)  of 

Mack,  H.  (1972).  Spatial  coherence  of  surface  waves.  Report  No.  SAAC-8,  Tele¬ 
dyne  Geotech,  Alexandria,  Virginia. 
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Figure  6.  Coherency  at  a  period  of  21.3  sec  (.0469  Hz)  versus  seismometer 
separation  for  the  Rayleigh-wave  recordings  of  Figure  5  (after 
Mack,  1972). 


Y2CT)  »  [1  +  cos  2Ttr  (kL  “  k2)]/2 

This  equation  would  fit  the  actual  data  of  Figure  6  nearly  as  well  as  the 
Gaussian  predictions  shown  there. 

In  order,  then,  to  define  the  typical  Rayleigh  signal  wavefield,  a  more 
detailed  analysis  of  the  observations  is  required.  Information  important  in 
determining  the  best  signal  model  is  contained  in  recording  phases  at  each 
array  sensor;  such  information  is  suppressed  in  the  simple  coherence  computa¬ 
tions,  but  it  can  be  recovered  with  alternative  methods  (discussed  later).  To 
overcome  NORSAR's  dimensions,  which  are  inadequate  to  fully  analyze  the  wave- 
field  nature,  signal  analysis  over  a  long  linear  array,  such  as  in  Figure  1, 
will  be  necessary.  Before  proceeding  to  the  data  analysis,  Chernov's  work, 
which  is  appropriate  to  a  wave  propagating  in  three  dimensions,  will  be  modified 
to  obtain  predictions  for  a  Rayleigh  wave,  which  essentially  propagates 
in  two  dimensions  near  the  earth's  surface. 
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CHERNOV’S  SCATTERING  THEORY  MODIFIED  FOR  RAYLEIGH  WAVES 


Chernov  (1962)  derived  analytic  expressions  that  statistically  described 
the  behavior  of  monochromatic,  acoustic,  plane  waves  as  they  travelled  through 
a  non-turbulent  fluid  medium  characterized  by  spatial  fluctuations  of  the  wave 
velocity  C  which  are  both  small  and  random.  Certain  simplifying  assumptions 
must  be  made  to  apply  these  results  to  the  problem  of  waves  propagating  across 
the  surface  of  a  solid,  such  as  seismic  Rayleigh  waves.  Complications  intro¬ 
duced  by  the  elliptical  particle  motion  of  Rayleigh  waves  are  ignored,  and  all 
displacements  are  assumed  to  be  along  the  direction  of  wave  propagation  in  the 
plane  of  the  earth's  surface.  This  assumption  not  only  has  the  effect  of 
replacing  Rayleigh  waves  by  acoustic  waves,  but  also  eliminates  the  considera¬ 
tion  of  effects  due  to  LR-LQ  mode  conversion  when  the  waves  are  scattered. 
Ignoring  the  dependence  of  the  waves  on  the  depth-varying  properties  of  the 
medium  below  the  surface  further  simplifies  this  problem.  The  problem  is  thus 
considered  strictly  two-dimensional,  and  Chernov's  theory  must  be  modified. 
Since  our  data  reduction  employs  spectral  analysis,  which  allows  analysis  of 
waves  within  only  a  narrow  band  of  frequencies,  the  assumption  that  the  waves 
are  monochromatic  is  acceptable.  Our  treatment  of  the  spatial  variation  of 
phase  velocity  C(x,y),  therefore,  does  not  need  to  involve  dependence  upon  the 
frequency  C  (x,y,w)  characteristic  of  dispersive  waves.  In  Chernov's  fluid 
medium  the  velocity  fluctuations  vary  in  time  as  well  as  in  space,  and  the  wave 
propagation  is  described  in  terms  of  averages  over  (an  infinite)  time  for  a 
given  source  and  a  given  receiver.  Although  in  the  solid  earth  the  velocity 
fluctuations  are  deterministic  and  permanent,  they  may  be  treated  in  Chernov's 
manner  if  averages  over  the  azimuths  of  incident  waves  from  (all  possible)  dif¬ 
ferent  sources  to  a  given  receiver  are  taken.  This  method  of  averaging  is  dis¬ 
cussed  more  fully  in  the  memorandum  of  Rivers  and  von  Seggern  (1978),  which  for 
convenience  is  reproduced  as  an  appendix  to  this  report. 


The  scattering  properties  of  an  inhomogeneous  medium  are  determined  by 
the  spatial  fluctuations  of  the  index  of  refraction,  which  are  defined  by 


p(x,y) 


C 

o 


C  (x,y) 


-1 


Rivers,  D.  W.  and  D.  H.  von  Seggern  (1978).  Theory  of  wave  propagation  in  a 
two-dimensional  random  medium,  Technical  Memorandum,  Teledyne  Geotech, 
Alexandria,  Virginia. 


where  Cq  is  the  average  phase  velocity  of  the  waves  in  the  medium.  Denoting 
the  mean  square  fluctuation  of  the  index  of  refraction  by  y^,  the  correlation 
coefficient  of  y  is  defined  to  be 

v  =  _ 

VJ(x1,y1)  p(x2,y2)  /  (6) 

If  the  fluctuation  y(x,y)  are  truly  random,  the  average  denoted  by  the  overbar 
in  (6)  depends  only  upon  the  distance  between  the  points  (x^,  y  )  and  (x2>  y2) 
and  not  upon  the  coordinates  of  the  points  themselves;  i.e.,  v=v(r).  In  the 
modification  of  Chernov's  theory  presented  here,  it  is  assumed  throughout  that 
the  form  of  the  correlation  coefficient  is  given  by  the  Gaussian  function 

v(r)  =  exp  (-r2/a2) 

where  the  "correlation  distance"  a.  represents  the  typical  linear  dimension  of 
the  inhomogeneities  which  give  rise  to  the  scattering.  Chernov's  theory  is  a 
statistical  method  of  describing  scattering  phenomena  within  a  medium  character¬ 
ized  by  the  parameters  y"^  and  a. 

Table  I  gives  theoretical  expressions  derived  in  Rivers  and  von  Seggern 
(1978)  for  seven  basic  quantities  which  characterize  the  scattering  process. 
These  seven  quantities  are: 

the  mean-square  fluctuation  of  the  amplitude  (in  natural  logarithmic 
units)  of  the  wave  at  a  given  receiver  relative  to  the  average  of 
the  amplitudes  at  all  receiver; 

the  mean-square  fluctuation  of  the  phase  (in  radians)  of  the  wave 
at  the  same  receiver  relative  to  the  phase  of  the  (hypothetical) 
average  wavefront; 

the  cross-correlation  coefficient  between  the  amplitude  and  phase 
fluctuations  at  the  same  receiver; 


B2  , 


P. 


R^(AL)  and  Rg  (AL) ,  the  spatial  autocorrelation  coefficients  for  amplitude 
or  phase  fluctuations  measured  at  two  different  receivers  whose 
separation  distance  AL  lies  along  the  direction  of  propagation  of 
the  incoming  wave;  and 

r^  (£)  and  Rg  (£)  the  corresponding  autocorrelation  coefficients  for  two 
receivers  whose  separation  £  lies  in  the  direction  transverse  to 
the  wave  propagation. 
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TABI.F.  t 


As  the  table  shows,  these  seven  quantities  can  be  expressed  as  functions  of 
the  four  variables:  p^",  a,  the  wavenumber  k,  and  L,  which  is  the  length  of 
of  wave's  path  through  the  medium.  For  a  narrow  spectral  band,  k  is  nearly 
constant.  Rivers  and  von  Seggern  (1978)  showed  that  variations  in  the  source 
to  receiver  distance  L  will  not  strongly  effect  our  analysis.  Therefore,  this 
study  will  determine  the  values  of  p7  and  a  for  those  inhomogeneities  in  the 
earth's  crust  causing  the  random  scattering  of  Rayleigh  waves  (near  20  sec 
period).  This  study  will  also  determine  the  applicability  of  Chernov’s  theory 
to  the  quantitative  prediction  of  fluctuations  due  to  that  scattering,  particu¬ 
larly  (AL)  and  R^  ( it)  . 

Note  the  conditions  that  govern  the  validity  of  the  result  in  Table  I. 

Tae  two-dimensional,  longitudinal-wave  model  of  Rayleigh  waves,  which  is  em¬ 
ployed  in  our  modification  of  Chernov's  theory,  has  already  been  described; 
also  the  assumptions  that  permit  replacement  of  equation  (6)  with  (7)  were 
stated.  One  other  set  of  Chernov's  key  assumptions  pertained  to  the  smallness 
of  fluctuations  in  density,  phase  velocity,  and  refractive  index: 


Ap<<Po>  Ac<<cq,  and  hence  p2<<l»  (8) 

The  validity  of  these  assumptions  must  be  tested  after  the  value  of  pr2"  is 
determined.  Further,  we  assume  that  the  scattering  caused  by  density  fluctu¬ 
ation  can  be  neglected  when  compared  with  scattering  caused  by  fluctuations 
in  phase  velocity.  This  assumption  is  valid  if,  for  a  wave  travelling  along 
the  x-axis. 


Zkp 


3(Ap) 

dX 


« 


Ac 


For  a  fairly  uniform  crustal  structure,  3(Ap)/3x  'v  01.01/gm/cc/km  would  be  a 
liberal  value  while  Ac/c^  'v  0.05  would  be  a  reasonable  value.  Using  p^  'v  3 
gm/cc  and  k  "v  0.1  for  20  sec  periods,  the  relation  is  satisfied. 

Another  assumption,  reflected  in  Table  I,  that  Chernov  used  was  that 
the  Born  approximation  was  valid  for  the  Rayleigh  waves.  This  first-order 
approximation,  which  Chernov  referred  to  as  the  method  of  small  perturbations, 
consists  of  substituting,  within  the  integrand  of  an  integral  equation,  the 
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expression  for  the  composite  wave  field  after  scattering  in  place  of  the  expres¬ 
sion  for  the  unscattered  incident  wave.  The  Born  approximation  is,  therefore, 
essentially  the  assumption  that  the  scattering  process  results  in  a  small 
perturbation,  i.e.,  that  the  amplitude  and  phase  fluctuations  are  both  small: 

«1  and  jA4>|«l. 

Ao' 

Chernov  claimed  that  his  formulas  were  valid  under  conditions  less  restrictive 
than  those  imposed  by  (10);  i.e.,  that  the  amplitude  and  phase  fluctuations  be 
small  over  the  distance  of  a  wavelength  A  rather  than  over  the  entire  source- 
to-receiver  path  length  L.  This  claim,  based  upon  Rytov's  method,  was  disputed 
by  Aki  (1973)  who  claimed  that  incident  wave  field  attenuation  due  to  energy 
lost  by  scattering  throughout  the  path  length  L  must  be  small.  Thus,  determin¬ 
ing  whether  the  Rayleigh-wave  data  is  consistent  with  the  assumption  of  small 
fluctuations  (10)  is  important. 

Chernov  also  assumed  that  the  source-to-receiver  path  length  L  is  large 
when  compared  to  the  correlation  distance,  which,  in  turn,  is  assumed  to  be 
large  when  compared  with  the  wavelength: 

L»a» 

k 

Chernov  defined  the  "wave  parameter"  as 


If  D  <<  1,  the  ray  equation  may  be  used  to  approximate  the  equations  of  motion; 
if  D  >>  1,  Fraunhofer  diffraction  theory  is  appropriate.  In  either  case  the 
condition  of  small  wavelength  (known  as  the  Fresnel  approximation),  ka  >>  1, 
must  be  obeyed.  Since  for  Rayleigh  waves  with  periods  of  20  sec  and  phase 
velocities  of  3.6  km/sec,  k  ,v»  .1  km  \  this  condition  implies  that  a^  should 
be  much  larger  than  10  km.  Because  this  condition  may  not  be  met,  there  are 
listed  in  Table  I,  wherever  appropriate,  two  expressions  for  the  given  func¬ 
tions.  The  more  precise  (and  more  cumbersome)  of  the  expressions  is  less 
sensitive  to  the  restriction  ka  >>  1  than  the  other.  Note,  however,  that  even 
the  more  exact  expressions  are  of  questionable  validity  where  the  size  of  the 
scatterers  is  approximately  equal  to  a  wavelength,  that  is,  kasl. 
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The  Fresnel  approximation  implies  that  most  scattering  deflects  the  waves 
forward  into  a  cone  of  aperture  angle  suggesting  that  the  direction  of 
the  scattered  wave  deviates  little  from  the  incident  wave.  Our  justification 
for  the  validity  of  this  approximation  is  based  upon  Capon's  (1970)  observa¬ 
tions  that  most  of  the  energy  in  a  typical  incident  wave  is  concentrated  in  a 
small  range  of  azimuths.  Finally,  note  that  even  if  ka  >>  1,  the  transition  to 
geometrical  acoustics  (k  -»•  °°)  is  not  allowed  by  the  Chernov  theory. 


RAYLEIGH-WAVE  SCATTERING  AT  NORSAR 


Data 

Measurements  were  made  of  the  amplitude  and  phase  of  20-sec  Rayleigh  waves 
recorded  at  NORSAR  for  a  selected  suite  of  ten  events  (listed  in  Table  II)  to 
provide  a  quantitative  description  of  fluctuation  magnitudes  and  of  the  corre¬ 
lation  among  the  fluctuation  at  different  locations  within  a  seismic  sensors 
array.  Rivers  and  von  Seggern  (1978)  explained  that  the  present  adaptation  of 
Chernov's  theory  utilizes  averages  over  source-to-receiver  paths  of  the  full 
360°  range  of  incident  azimuths.  However,  this  procedure  could  not  be  fully 
implemented  because  of  the  location  of  suitable  seismic  regions.  Still,  the 
events  chosen  for  analysis  in  this  study  were  selected  to  attain  as  wide  a 
range  as  possible  and  yield  even  spacing  in  incident  azimuths.  Nevertheless, 
the  range  of  (geometrical)  azimuths,  shown  in  Table  II,  is  only  about  135°. 

The  range  is  limited  because  of  the  need  to  include  only  events  for  which  the 
great-circle  paths  to  NORSAR  did  not  cross  any  major  source  of  non-random  scat¬ 
tering,  such  as  an  ocean-continent  boundary.  Also,  selection  was  restricted 
because  we  wanted  events  that  took  place  within  as  short  a  time  period  as  pos¬ 
sible  to  minimize  the  effects  of  potential  long-term  instabilities  in  the 
response  of  the  instruments  at  NORSAR.  To  check  measurement  repeatability,  one 
event  was  chosen  that  was  an  aftershock  of  another;  thus,  two  events  in  the 
chosen  data  set  occurred  at  approximately  the  same  epicenter  within  three  days 
of  each  other.  Finally,  an  effort  was  made  to  pick  events  with  as  large  a  Ray¬ 
leigh-wave  magnitude  as  possible  so  that  most  of  the  measured  fluctuations 
would  stem  from  wave  scattering  and  not  from  background  noise. 

Computational  Methods 

The  data  for  each  event  consisted  of  seismograms  from  twenty-two  long- 
period  vertical  instruments  within  the  array.  These  twenty-two  seismograms 
were  analyzed  in  frequency-wavenumber  space  using  the  program  that  Smart  (1972) 
developed,  indicating  the  time  window  (chosen  to  be  256  sec  long)  where  the 
signal  power  was  greatest,  the  phase  velocity,  and  the  incident  azimuth,  all 
as  a  function  of  frequency.  The  FKCOMB  results  confirmed  that  no  strong 


Smart,  E.  (1972).  FKCOMB,  a  fast  general-purpose  array  processor.  Report  No 
SAAC-9,  Teledyne  Geotech,  Alexandria,  Virginia. 
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multipath  arrivals  contaminated  the  data  for  the  selected  events.  The  time 
windows  chosen  with  the  aid  of  FKCOMB  as  the  best  for  waves  of  19.7  sec  peri¬ 
od  were  then  used  to  compute  the  amplitude  and  phase  spectra  for  each  of  the 
twenty-two  seismograms;  for  any  of  the  events,  this  time  window  was  identical 
for  all  22  sensors.  We  then  determined  the  amplitude  at  each  instrument  by 
smoothing  the  natural  logarithms  of  the  spectral  amplitudes  for  the  periods 
21.3,  19.7,  and  18.3  sec,  using  Hanning  weights  of  1/4,  1/2,  1/4,  respectively. 
The  phase  of  the  wave  at  each  intrument  was  assumed  to  be  the  spectral  phase 
for  the  period  19.7  sec. 

Amplitude  fluctuations  for  each  event  at  each  instrument  were  then  calcu¬ 
lated  from  the  values  found  for  each  instrument.  Also  for  each  event  a  grid 
search  of  azimuths  and  a  linear  least-squares  fit  to  the  phase  yielded  the 
azimuth  and  phase  velocity  of  the  plane  wave  best  fitting  the  spectral  phase 
data.  The  azimuths  and  phase  velocities  of  the  best-fitting  wave  found  with 
this  procedure  were  consistent  with  values  calculated  by  FKCOMB.  These  phase 
fluctuations  were  assumed  to  be  the  advance  or  the  delay  of  the  spectral  phase 
at  each  instrument  relative  to  the  phase  of  the  hypothetical  plane  wave  as  it 


passes  through  that  point 

in  the 

array.  Thus, 

we  have 

B.  .  =  AnA. . 

inA. 

i 

i  =  1,2,  .. 

.,  10 

j  -  1,2, 

...,  22 

(13) 

S.  .  =  <J>.  (j)  - 
lj  1  J 

i  =  1,2,  .. 

.,  10 

j  =  1,2, 

...,  22 

a*) 

where  denotes  the  phases  of  the  (hypothetical)  coherent  incident  wave  as 

measured  at  the  jth  instrument.  Note  that  <(>^(j)  is,  of  course,  dependent  upon 
values  assumed  for  the  phase  velocity  and  azimuth  of  the  incident  wave.  The 
phase  velocities  computed  for  the  best-fitting  plane  waves  were  not  the  same 
for  each  event;  in  fact,  they  varied  from  3.503  to  3.976  km/sec  (see  Figure  7). 
Since  theoretical  models  of  the  crustal  structure  underneath  NORSAR  (Massd 
and  Alexander,  1974)  yield  a  phase  velocity  (for  waves  of  19.7  sec  period)  of 
3.64  km/sec  and  since  this  value  should  presumably  be  constant  for  all  events, 
we  decided  to  proceed  with  the  analysis  in  two  different  ways. 


Masse,  R»  P*,  and  S.  S.  Alexander  (1974).  Compressional  velocity  distribution 
beneath  Scandinavia  and  western  Russia,  Geophys.  J. ,  39 ,  587-602. 
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One  mode  of  analysis  calculated  on  the  basis  of  a  constant  phase  velo¬ 

city  of  3.64  km/sec  for  each  event;  the  other  mode  used  the  phase  velocity 
which  best  fit  the  data  for  the  ith  event.  Our  data  base  thus  consisted  of 
the  set  of  amplitude  residuals  B^.  and  two  sets  of  phase  residuals  . 


Results 

These  three  sets  of  residuals  are  listed  in  Tables  III,  IV,  and  V.  In 
Table  III,  for  the  amplitude  fluctuations,  the  extreme  right-hand  column  con¬ 
tains  for  each  instrument  the  mean  fluctuation  and  its  standard  deviation  over 
all  events.  (Note  that  not  every  instrument  was  operative  during  each  of  the 
ten  events.)  For  every  instrument  the  95 %  confidence  limits  (±2  standard 
deviation/number  of  events)  on  the  mean  overlaps  zero  demonstrating  that  no 
significant  net  amplitude  biases  existed  that  had  to  be  removed  from  the  data 
for  each  instrument.  At  the  bottom  of  each  event  column  in  Table  III  the 
standard  deviation  of  the  amplitude  fluctuations  for  that  event  over  all  the 
instruments  is  listed.  (The  mean  of  the  fluctuations  is  zero  by  definition 

-7T  *  1^ 

in  equation  (13).)  This  value  is  equal  to  (B^  )  2  for  the  ith  event.  The  aver- 

i 

age  of  this  quantity  for  all  ten  events  is  0.249,  corresponding  to  a  standard 
deviation  in  surface-wave  magnitude  of  (0.249)/  jn  (10)  =  0.108.  A  typical 

T  2 

value  of  B  is  thus  (0.249)  =  0.062.  The  rms  fluctuation  of  0.249  is  only 

in  fair  agreement  with  our  hypothesis  of  small  fluctuation  in  amplitude.  Sim¬ 
ilar  results  hold  for  Tables  IV  and  V,  which  list  the  two  sets  of  phase  resi¬ 
duals  S...  The  mean  fluctuation  over  all  events  and  all  instruments  had  a 
standard  deviation  of  about  13  degrees  for  both  data  sets,  which  corresponds 

— *2  ^  j 

to  a  value  of  S  =  0.051  in  radian  measure.  The  rms  (S^)'2  =  0.227,  which 
again  is  somewhat  larger  than  the  fluctuations  considered  in  Chernov's  theory. 


The  repeatability,  and  to  a  certain  extent  the  reliability,  of  the  mea¬ 
surements  can  be  determined  by  comparing  the  fluctuations  for  the  two  events 
in  the  same  aftershock  sequence,  events  600  and  700.  Since  the  epicenter-to- 
NORSAR  path  for  these  two  events  is  nearly  identical,  the  residuals  at  each 
instrument  should  also  be  identical,  an  observation  supported  in  an  examina¬ 
tion  of  Tables  III,  IV,  and  V.  (One  glaring  exception,  however,  is  instrument 
number  16,  whose  amplitude  residual  for  event  600  suggests  a  malfunction.) 
Although  most  other  residuals  for  the  two  events  are  nearly  the  same,  they  are 
not  exactly  so.  This  difference  is  thought  to  stem  from  slight  differences  in 
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Amplitude  Fluctuations  at  NORSAR 
(natural  log  of  amplitudes,  with  event  mean  removed) 
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*  ***  -rM 


Phase  Fluctuations  at  NORSAR  Assuming  Phase  Velocity  *  3.641  km/sec 
^fluctuation  *  calculated  -  observed  residual  in  degrees) 


+  method  of  computing  fluctuations  for  a  given  event  requires  that  mean  =0.0 
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absolute  value  of  the  fluctuation  at  this  seismometer  is  greater  than  twice  the  standard  deviation 
for  this  event 


the  paths  of  the  incident  waves  and  not  from  variations  within  a  three-day 
period  of  the  instrument  responses.  If  this  were  the  case  it  would  strongly 
suggest  that  instrument  responses,  at  least  for  some  instruments,  were  subject 
to  significant  random  fluctuations  during  the  period  the  data  set  represents. 


Rivers  and  von  Seggern  (1978)  explained  that  important  parameters  in  the 

2  2 

Chernov  theory  are  the  ratio  of  B  to  S  and  the  correlation  coefficient  p 
between  the  amplitude  and  phase  fluctuations.  Capon  (1974)  derived  expressions 
for  the  confidence  limits  about  the  theoretical  values  of  B^/S^- and  p  within 
which  estimates  of  those  quantities  should  fall  based  upon  N  observations. 

The  confidence  limits  R  about  p  are  given  by 

R  =  tanh  (Z  +  2a)  (15) 

where  . 

Z  =  1/2  in  (  J  I  g  )  (16) 

and 

o  =  (N  -  3)  "1/2  (17) 


We  surround  our  observed  values  of  p  with  error  bars  by  substituting  the 
observed  value  in  place  of  the  theoretical  value  in  (16).  Error  bars  about 
2  2 

B  /s  are  obtained  by  noting  that  this  quantity  has  the  F  distribution.  The 

2  2 

values  of  B  /S  ,  p>  and  their  95%  confidence  limits  are  given  in  Table  VI. 

2  2 

For  every  event  the  ratio  B  /S  is  greater  for  the  data  set  where  the  phase 
velocity  of  the  best-fitting  incident  wave  rather  than  the  constant  C  =  3.64 
km/sec  was  used,  because  the  best-fitting  phase  velocity  is  defined  to  be  the 
one  for  which  S^"  is  a  minimum  for  that  event.  The  difference  between  the 
values  of  p  found  by  using  the  two  different  sets  of  phase  residuals  is  impos¬ 
sible  to  predict. 

TF 

The  phase  residuals  and,  therefore,  the  quantities  B  /S  and  p  are  depen¬ 
dent  not  only  upon  the  assumed  phase  velocity  of  the  incident  wave  but  also 
upon  the  value  assumed  for  its  azimuth.  Table  II  shows  that  the  azimuths 
found  by  FKCOMB,  which  minimized  the  residuals  s"^  f°r  a  given  event,  did  not 
always  agree  with  the  geometrical  NORSAR-to-epicenter  azimuth  for  that  event. 
The  sensitivity  of  and  p  to  variations  in  the  assumed  azimuth  of  the  inci¬ 
dent  wave  should  be  examined.  This  is  done  for  event  200  in  Table  VII,  which 
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TABLE  VI 


Ratios  and  correlations  of  single-receiver 
amplitude  and  phase  fluctuations  at  NORSAR 

Complete  Data  Set 

Phase  ^2 

Azimuth  velocity  -  Confidence  Confidence 


Event 

(degrees) 

(km/sec) 

S2 

interval 

P 

interval 

200 

106.9 

3.641 

0.448 

0.185 

to  1.084 

-0.385 

-0.700 

to  0.056 

3.550 

0.701 

0.290 

1.696 

-0.081 

-0.495 

0.363 

300 

106.7 

3.641 

0.591 

0.240 

1.454 

0.227 

-0.240 

0.608 

3.550 

0.641 

0.261 

1.577 

0.246 

-0.220 

0.621 

400 

80.9 

3.641 

0.280 

0.114 

0.689 

0.616 

0.238 

0.832 

3.829 

0.355 

0.144 

0.873 

0.290 

-0.174 

0.650 

500 

16.0 

3.641 

0.490 

0.194 

1.240 

-0.339 

-0.687 

0.136 

3.575 

0.512 

0.202 

1.295 

-0.306 

-0.668 

0.172 

600 

73.0 

3.641 

2.473 

1.005 

6.084 

0.323 

-0.139 

0.670 

3.720 

2.835 

1.152 

6.974 

0.267 

-0.199 

0.634 

700 

73.9 

3.641 

1.702 

0.692 

4.187 

0.538 

0.125 

0.792 

3.707 

1.832 

0.745 

4.507 

0.492 

0.064 

0.768 

800 

62.8 

3.641 

1.156 

0.484 

2.763 

-0.283 

-0.630 

0.157 

3.668 

1.196 

0.500 

2.653 

-0.348 

-0.671 

0.086 

900 

152.4 

3.641 

0.507 

0.212 

1.212 

0.349 

-0.085 

0.672 

3.579 

0.519 

0.217 

1.240 

0.312 

-0.126 

0.648 

1000 

6.0 

3.641 

1.405 

0.571 

3.456 

-0.192 

-0.585 

0.274 

3.503 

1.510 

0.614 

3.715 

-0.135 

-0.545 

0.327 

1100 

85.4 

3.641 

1.432 

0.582 

3.523 

-0.100 

-0.519 

0.358 

3.976 

4.262 

1.733 

10.485 

-0.529 

-0.787 

-0.113 
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TABLE  VII 


Effect  of  varying  the  assumed  NORSAR-to-epicenter  azimuth  for  event  200 
Phase  velocity  assumed  to  be  3.641  km/sec. 


Complete  Data  Set 


Assumed  azimuth 
(degrees) 


20  =  sum  of  squares  of 

phase  residuals  (degrees^) 


96.9 

14,740 

0.025 

-0.259 

97.9 

12,030 

0.031 

-0.264 

98.9 

9,611 

0.038 

-0.270 

99.9 

7,482 

0.049 

-0.280 

100.9 

5,644 

0.065 

-0.293 

101.9 

4,099 

0.090 

-0.311 

102.9 

2,849 

0.129 

-0.336 

103.9 

1,895 

0.194 

-0.368 

104.9 

1,238 

0.297 

-0.405 

105.9 

880 

0.418 

-0.424 

106.9 

821 

0.448 

-0.385 

107.9 

1,061 

0.347 

-0.294 

108.9 

1,603 

0.230 

-0.206 

109.9 

2,446 

0.150 

-0.142 

110.9 

3,591 

0.102 

-0.098 

111.9 

5,038 

0.073 

-0.069 

112.9 

6,788 

0.054 

-0.048 

113.9 

8,841 

0.042 

-0.034 

114.9 

11,200 

0.033 

-0.024 

115.9 

13,860 

0.027 

-0.016 

116.9 

16,820 

0.022 

-0.011 
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assumes  the  phase  velocity  of  the  incident  wave  to  be  3.64  km/sec  and  in  Table 

VIII,  which  uses  the  best-fitting  phase  velocity  for  each  value  assumed  for 

_  _ - 

the  incident  azimuth.  Clearly  and  B  /S^  change  markedly  for  azimuths  only 

2 

a  few  degrees  from  the  values  for  which  S  is  a  minimum.  The  correlation  p  also 
exhibits  a  strong,  but  less  dramatic,  dependence  upon  azimuth.  The  values 
~2~  2~ 

B  /S  and  p  for  event  200  would  have  been  significantly  different  if  the  geo¬ 
metrical  azimuth  (99.0°)  rather  than  the  value  found  by  FKCOMB  (106.9°)  were 
used. 

Table  III,  IV,  and  V  show  several  amplitude  and  phase  fluctuations  whose 
magnitude  exceeds  twice  the  standard  deviation  of  fluctuations  for  that  event. 
That  is,  for  certain  events  i  and  certain  instruments  j, 

|Bijj>2(B.1)  1/2  (18) 

and  similarly  for  S...  Because  these  large  (and  possibly  spurious)  fluctuations 

2  — r 

have  a  signicant  effect  upon  B  ,  S  >  and  P*  another  data  set,  the  "edited" 
set,  was  formed  by  omitting  these  anomalous  values  from  the  complete  data  set. 

As  already  noted,  we  also  wanted  the  data  set  to  reflect  as  wide  and  as  well¬ 
spaced  a  range  of  azimuths  as  possible  so  the  statistical  average  over  source- 
to-receiver  paths  in  the  Chernov  theory  would  be  valid;  for  this  reason  one  of 
the  events  in  the  aftershock  sequence  (event  600)  was  not  included  in  the 
edited  data  set.  Ths  omission  prevents  excess  emphasis  on  a  particular  path 
in  the  averaging  process.  (Athough  event  200  and  300  had  almost  the  same 
FKCOMB  azimuths  for  the  incident  wave,  their  geometrical  NORSAR-to-epicenter 
azimuths  are  different.) 

For  the  edited  data  set  (Tables  IX,  X,  and  XI) ,  the  rms  amplitude  fluc¬ 
tuation  is  0.147  natural  log  unit,  a  figure  equivalent  to  a  standard  devia¬ 
tion  of  0.064  in  the  measurement  of  the  Rayleigh-wave  magnitude  Mg  from  20 
sec  waves.  The  rms  phase  fluctuation,  calculated  by  assuming  a  phase  velocity 
of  3.64  km/sec,  is  11.2  degrees,  or  0.195  radians;  the  fluctuation  calculated 
on  the  basis  of  the  best-fitting  phase  velocity  is  9.7  degrees,  or  0.169  ra¬ 
dians.  The  amplitude  fluctuations,  therefore,  are  much  smaller,  and  the  phase 
fluctuations  somewhat  smaller,  for  the  edited  data  set  than  for  the  complete 
data  set.  The  values  of  B and  P  calculated  for  the  edited  data  set  are 
listed  in  Table  XII. 
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TABLE  VIII 


Effect  of  varying  the  assumed  NORSAR-to-epicenter  azimuth  for  event  200 
Best-fitting  phase  velocity  used  for  each  azimuth 


Complete  Data  Set 


Assumed 

azimuth 

(degrees) 


o 

20  S  =  sum  of  squares 
of  phase  residuals 
_ (degrees^) _ 


phase 
velocity 
(km/ sec) 


96.9 

14,740 

0.025 

-0.261 

3.644 

97.9 

12,030 

0.031 

-0.254 

3.629 

98.9 

9,589 

0.038 

-0.246 

3.616 

99.9 

7,432 

0.050 

-0.239 

3.603 

100.9 

5,559 

0.066  „ 

-0.232 

3.592 

101.9 

3,976 

0.093 

-0.225 

3.582 

102.9 

2,686 

0.137 

-0.217 

3.573 

103.9 

1.693 

0.217 

-0.208 

3.565 

104.9 

1,000 

0.368 

-0.193 

3.559 

105.9 

610 

0.604 

-0.156 

3.554 

106.9 

525 

0.701 

-0.081 

3.550 

107.9 

748 

0.492 

-0.003 

3.547 

108.9 

1,280 

0.287 

0.040 

3.546 

109.9 

2,122 

0.173 

0.059 

3.545 

110.9 

3,275 

0.112 

0.066 

3.546 

111.9 

4,739 

0.078 

0.066 

3.549 

112.9 

6,513 

0.056 

0.063 

3.552 

113.9 

8,597 

0.043 

0.058 

3.557 

114.9 

10,990 

0.033 

0.052 

3.564 

115.9 

13,690 

0.029 

0.044 

3.571 

116.9 

16,690 

0.022 

0.037 

3.580 
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Amplitude  fluctuations  at  NORSAR  (natural  log  of  amplitudes,  with  event  mean  removed) 
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Average  fluctuation  for  all  seismometers  and  all  event: 
0.00  +  .147 


Phase  fluctuations  at  NORSAR  assuming  phase  velocity  =  3.641  km/sec 
(fluctuation  =  calculated-observed  residual  in  degrees) 
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Average  flucutuation  for  all  seismometers  and  all  event: 


TABLE  XII 


Ratios  and  correlations  of  single-receiver 
amplitude  and  phase  fluctuations  at  NORSAR 


Edited  Data  Set 


Event 

Azimuth 

(degrees) 

Phase 

velocity 

(km/sec) 

B2 

S2 

Confidence 

interval 

P 

Confidence 

interval 

200 

106.9 

3.641 

0.425 

0.168 

to  1.075 

-0.359 

-0.699 

to  0.113 

3.548 

0.715 

0.283 

1.809 

-0.096 

-0.527 

0.375 

300 

106.7 

3.641 

0.437 

0.173 

1.106 

-0.038 

-0.484 

0.424 

3.559 

0.473 

0.187 

1.197 

-0.011 

-0.463 

0.445 

400 

80.9 

3.641 

0.465 

0.184 

1.176 

0.530 

0.100 

0.793 

3.772 

0.576 

0.228 

1.457 

0.213 

-0.267 

0.608 

500 

16.0 

3.641 

0.490 

0.194 

1.240 

-0.339 

-0.687 

0.136 

3.575 

0.512 

0.202 

1.295 

-0.306 

-0.668 

0.172 

700 

73.9 

3.641 

1.754 

0.675 

4.560 

0.196 

-0.298 

0.608 

3.686 

1.864 

0.717 

4.846 

0.144 

-0.346 

0.573 

800 

62.8 

3.641 

0.695 

0.287 

1.682 

-0.003 

-0.434 

0.429 

3.687 

0.776 

0.321 

1.878 

-0.075 

-0.491 

0.369 

900 

152.4 

3.641 

0.689 

0.272 

1.743 

0.499 

0.058 

0.777 

3.587 

0.707 

0.279 

1.789 

0.477 

0.029 

0.765 

1000 

6.0 

3.641 

1.405 

0.571 

3.456 

-0.192 

-0.585 

0.274 

3.503 

1.510 

0.614 

3.715 

-0.135 

-0.545 

0.327 

1100 

85.4 

3.641 

0.355 

0.137 

0.923 

0.379 

-0.106 

0.719 

4.000 

2.705 

1.040 

7.033 

0.180 

-0.313 

0.597 
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While  the  logarithmic  amplitude  fluctuation  mentioned  above  was  calculated 
by  smoothing  the  spectral  amplitudes  for  the  three  distinct  periods  21.3,  19.7, 
and  18.3  sec,  the  phase  fluctuation  was  calculated  upon  the  basis  of  spectral 
phases  for  the  distinct  period  19.7  sec  alone.  A  third  data  set  was  created 
that  treated  each  of  these  three  periods  separately  by  calculating  unsmoothed 
logarithmic  amplitudes  and  phase  delays  on  the  basis  of  each  period's  data 
alone.  As  before,  two  sets  of  phase  residuals  were  calculated;  for  one  set 
the  phase  velocities  were  assumed  constant  for  all  events.  For  the  other  set, 
the  phase  velocity  used  for  each  event  was  the  best-fitting  value  for  that 
event.  However,  it  was  assumed  that  the  incident  azimuth  was  the  same  for 
each  period,  namely  the  value  which  was  found  for  the  waves  of  period  19.7  sec. 
Only  events  and  instruments  which  were  included  in  the  "edited"  data  set  were 
included  in  the  third  data  set.  The  amplitude  and  phase  fluctuation  for  this 
data  set  are  typified  by  event  200,  shown  in  Table  XIII.  The  phase  residuals 
are  smaller  for  the  19.7  sec  waves  than  for  the  other  two  periods,  since  the 
chosen  azimuth  minimizes  the  residuals  for  this  period.  A  fourth  data  set 
was  created  like  the  third  except  for  each  period  the  azimuth  of  the  incident 
wave  was  assumed  to  be  the  value  that  minimized  the  phase  residuals  for  that 
period  alone.  The  fluctuations  in  this  fourth  data  set  (not  tabulated)  were 
the  same  size  as  those  in  the  third,  save  that  the  phase  residuals  for  periods 
of  21.3  and  18.3  sec  were  now  about  the  same  size  as  those  for  19.7  sec.  The 
third  data  set  thus  represents  the  same  assumption  as  the  "complete"  and  "edited" 
sets,  i.e.,  that  the  "true"  azimuth  of  the  waves  within  the  spectral  band  under 
consideration  is  given  by  the  value  that  minimizes  the  phase  residuals  for  the 
period  in  the  middle  of  that  band,  an  assumption  not  made  in  calculating  the 
fourth  data  set.  The  best-fitting  azimuths  and  phase  velocities  (Table  XIV 
and  Figure  8)  vary  significantly  within  the  chosen  spectral  band  for  several 
events. 

The  quantities  B^/s^  and  p  were  not  evaluated  for  each  event  and  each 

frequency  using  the  third  and  fourth  data  sets.  Instead,  these  data  sets, 

as  well  as  the  first  two,  are  now  used  in  a  different  form  of  analysis. 

Chernov's  theory  involves  taking  averages  over  all  possible  incident  azimuths; 

so,  average  values  of  B^»  S  ,  and  p  were  found  for  each  of  the  four  data  sets 

—j 

using  the  combined  data  from  all  events.  Thus,  instead  of  ten  values  of  B 
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Amplitude  and  phase  fluctuations  at  NORSAR  for  event  200 
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Azimuths  and  phase  velocities  yielding  best 
fit  to  data  for  each  of  three  frequencies 


Period  = 

21.3  sec 

Period  = 

19.7  sec 

Period  = 

18.3  sec 

Event 

0  (degrees) 

c (km/sec) 

0  (degrees) 

c(km/sec) 

0  (degrees) 

c (km/sec) 

200 

102.0 

3.643 

107.0 

3.543 

102.4 

3.601 

300 

129.4 

3.649 

106.4 

3.561 

110.3 

3.603 

400 

83.9 

3.618 

80.5 

3.770 

83.6 

3.431 

500 

20.4 

3.752 

15.5 

3.576 

14.0 

3.487 

700 

71.9 

3.631 

75.0 

3.679 

71.2 

3.553 

800 

66.0 

3.670 

62.5 

3.693 

61.4 

3.403 

900 

155.5 

3.489 

152.9 

3.593 

170.8 

3.543 

1000 

353.0 

3.337 

9.1 

3.504 

18.4 

3.369 

1100 

78.7 

3.751 

84.1 

4.012 

83.2 

3.658 
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being  computed  on  the  basis  of  as  many  as  twenty-two  observations  each,  one 
value  is  computed  on  the  basis  of  204  observations.  (Raw  data  from  diiferent 
events  must  not  be  used  together.  Rather,  use  the  fluctuations,  i.e.,  data 
from  which  had  been  removed  the  mean  of  the  amplitude  and  the  phase  delays  of 
the  incident  wave  for  each  event.) 

Table  II  indicates  that  the  azimuth  determined  by  FKCOMB  differs  marked! v 

from  the  geometrical  NORSAR-to-epicenter  azimuth  for  events  500  and  1000. 

Their  geometrical  azimuths  and  source  regions  suggest  that  the  eastern  and, 

especially,  the  northern  edges  of  the  Asian  continent  scatter  rays  from  these 

events.  Possibly  as  a  result  of  this  scattering,  the  amplitude  and  phase 

residuals  for  event  1000  are  significantly  greater  than  those  for  the  other 

events.  To  avoid  the  possible  systematic  (non-random)  scattering  phenomena 

~ 2  ~ 2 

that  created  these  fluctuations,  B  >  s  >  and  p  have  been  redetermined  from 
the  second  and  third  data  sets  by  deleting  the  data  from  these  two  events. 
Therefore,  six  data  sets  were  used  to  determine  the  average  over  several 
events  of  B  ,  S  ^  and  p.  addition,  each  of  these  data  sets  has,  in  turn, 
two  sets  of  phase  residuals,  one  set  assuming  the  constant  phase  velocity 
and  the  other  set  using  a  different  value  of  C  for  each  event.  These  results 
are  summarized  in  Table  XV. 

Wave  Parameter  D 

The  wave  parameter  D  is  an  important  quantity  in  the  Chernov  theory. 

Rivers  and  von  Seggern  (1978)  showed  that  D  can  be  determined  from  either 

/s^  or  from  P  by  means  of  the  formula  listed  in  Table  I.  Also,  the  wave 

parameter  can  be  eliminated  between  the  expressions  for  Bz/Sz  and  for  p,  thereby 

determining  one  of  these  quantities  as  a  function  of  the  other.  Following  Aki's 

(1973)  approach,  D  is  calculated  from  B7/?7  and  also  independently  from  p.  A 

value  which  lies  between  these  two  independent  values  is  found  by  determining 

the  point  on  the  theoretical  B^ /s^  versus  p  curve  (shown  in  Rivers  and  von  Seggern 

(1978))  closest  to  the  observed  point  (p,  B^/s^)>  this  value  is  designated 

D  .  Tables  XVI  and  XVII  list  these  three  values  of  D  for  each  event  for 
curve 

the  complete  and  edited  data  sets.  Another  value  of  D  was  also  calculated 

using  these  three  methods  for  each  of  the  six  data  sets  that  are  composites 

of  all  events;  it  was  from  composite  data  that  Capon  (1974)  found  his  value 

of  D  .  These  results  are  shown  in  Table  XVIII. 
curve 
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TABLE  XV 


Ratios  and  correlations  of  single-receiver 
amplitude  and  phase  fluctuations  at  NORSAR 


Average  over  all  events 


Data  Set 

Phase 

Velocity 

(km/sec) 

B2 

7 

P 

Number  of 

Detections 

Complete 

3.641 

1.028 

0.039 

204 

variable 

1.211 

0.006 

Edited 

3.641 

0.894 

-0.009 

172 

variable 

1.100 

-0.020 

Edited  , 

3.641 

0.586 

0.267 

133 

7  events 

variable 

0.838 

0.196 

Edited, 

3.705, 

3  separate 

3.641, 

frequencies 

3.587 

0.646 

0.010 

172 

variable 

0.683 

0.025 

Edited, 

3.705, 

3  separate 

3.641, 

frequencies. 

3.587 

0.569 

0.107 

133 

7  events 

variable 

0.608 

0.104 

Edited, 

3.705, 

3  separate 

3.641, 

frequencies , 

3.587 

1.353 

0.009 

172 

different 

azimuth  for 

variable 

1.536 

0.060 

each 

frequency 
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"-Sir.  W 


Values  of 


Event 

Phase 
velocity 
(km/  sec) 

200 

3.641 

3,550 

300 

3.641 

3.550 

400 

3.641 

3.829 

500 

3.641 

3.575 

600 

3.641 

3.720 

700 

3.641 

3.707 

800 

3.641 

3.668 

900 

3.641 

3.579 

1000 

3.641 

3.503 

1100 


3.641 


49 


TABLE  XVII 

Values  of  the  wave  parameter  D 
Edited  Data  Set 


Phase 

D 

7/s2 

Event 

velocity 

(km/sec) 

200 

3.641 

11.2 

3.548 

76 

300 

3.641 

12.0 

3.559 

14.4 

400 

3.641 

14.0 

3.772 

26 

500 

3.641 

16 

3.575 

18 

700 

3.641 

3.686 

800 

3.641 

68 

3.687 

130 

900 

3.641 

63 

3.587 

71 

1000 

3.641 

3.503 

1100 

3.641 

4.000 

7.8 

TABLE  XVIII 


Values  of  the  wave  parameter  D 
Average  over  all  events 


Data  Set 

Phase 

velocity 

(km/sec) 

Complete 

3.641 

variable 

Edited 

3.641 

variable 

Edited, 

3.641 

7  events 

variable 

Edited, 

3.705, 

3  separate 

3.641, 

frequencies 

3.587 

variable 

Edited, 

3.705, 

3  separate 

3.641, 

frequencies. 

3.587 

7  events 

variable 

Edited, 

3.705, 

3  separate 

3.641, 

frequencies , 
different 

3.587 

azimuth  for 
each 

frequency 

variable 

D 

— r  D  D 

B2/g2  p  curve 


830 


28  14  25 

310  32  140 


44 

61 


25  150  34 

32  160  43 


700 
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Discussion  of  Experimental  Observation 

Note  that  about  half  the  data  consists  of  values  with  D  undefined  because 

2  2 

the  observed  values  of  p,  or  B  /S  ,  lie  outside  the  range  of  permissable  values 
useful  in  calculating  D.  The  Chernov  theory  predicts  that 

0  <  P  <  0.5  (19) 


but  frequently  our  data  did  not  meet  these  conditions.  Perhaps,  the  inequality 
(20),  and  also  indirectly  (19),  should  not  always  be  expected  to  hold  because 
for  several  of  the  data  sets  efforts  were  made  to  minimize  the  phase  residuals 
by  selecting  the  best  value  of  the  incident  azimuth  and  phase  velocity.  This 
procedure  indicates  the  need  to  determine  the  proper  average  values  to  use 
when  calculating  amplitude  and  phase  residuals.  (So  far  we  have  assumed  that 
these  values  were:  for  amplitude,  the  average  of  the  logarithmic  amplitudes 
for  all  the  instruments;  and,  for  phase,  the  plane  wave  with  azimuth  and  vel¬ 
ocity  that  best  fits  the  phase  data.)  The  values  calculated  for  the  fluctua¬ 
tions  would  have  been  different  if  these  mean  values  were  chosen  differently. 

One  reason  that  they  perhaps  should  have  been  different  is  that  the  mean  values 
themselves  may  reflect  the  random  scattering  process.  For  example,  if  the  data 
set  had  consisted  of  magnitudes  measured  not  only  at  NORSAR,  but  also  at  other 
stations  nearby  the  array,  the  average  of  the  amplitudes  probably  would  have 
been  different  from  the  value  determined  by  NORSAR  data  alone.  Because  we 
think  that  some  of  this  difference  stems  from  random  scattering,  the  question 
arises  of  the  applicability  of  Chernov's  theory  to  random  fluctuations  about 
a  value  which  is  itself  a  random  fluctuation.  Again,  note  that  the  fluctuations 
are  assumed  small,  but  the  difference  between  the  average  magnitude  at  NORSAR 
and  the  average  magnitude  at  a  point  several  hundred  kilometers  away  can  be 
considerably  greater  than  the  average  magnitude  difference  within  the  hundred 
kilometer-wide  array.  The  fluctuations  measured  across  NORSAR,  then,  may 
represent  scattering  on  a  smaller  scale  than  that  for  which  Chernov's  theory 
is  valid,  given  the  true  size  of  the  inhoraogeneities  which  are  the  most  effi¬ 
cient  scatters  of  waves  of  70  km  wavelength.  These  points  will  be  further 
discussed  following  analysis  of  data  from  the  linear  array  in  the  Southwestern 
United  States. 
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Another  problem  in  analyzing  the  data  is  that  the  azimuth  and  phase  velo¬ 
city  for  the  best-fitting  plane  wave  can  vary  significantly  within  a  narrow 
spectral  band,  further  demonstrating  the  difficulty  of  choosing  the  mean 
value  about  which  fluctuations  should  be  measured.  This  variation  of  azimuth 
with  frequency  suggests  that  Rayleigh-wave  scattering  may  differ  from  the 
model  that  Aki  (1973)  and  Capon  (1974)  assumed  for  P-wave  scattering.  In 
their  model,  amplitude  and  phase  fluctuations  measured  at  LASA  stemmed  from 
scattering  within  the  crust  of  a  single  coherent  wave  incident  from  the  mantle. 
However,  because  of  the  Rayleigh  waves'  much  longer  path  through  the  scattering 
medium,  the  fluctuation  measured  at  NORSAR  results  from  scattering  through  the 
array  of  several  different  rays  incident  at  different  angles,  and  these  rays 
themselves  might  have  been  scattered  from  other  rays  rather  than  from  the  same 
wavefront.  If  this  multiple  scattering  occurs,  then  the  Born  approximation, 
which  separates  the  wave  field  into  an  incident  wave  carrying  most  of  the 
energy  and  a  diffracted  wave  of  small  amplitude,  breaks  down.  In  any  case, 
our  data  is  frequently  inconsistent  with  inequalities  (19)  and  (20)  predicted 
by  the  Chernov  theory;  however,  note  that  values  measured  for  B^/^  an^  f°r 
p  lie  within  large  confidence  limits. 

Spatial  Correlations 

The  spatial  correlation  among  fluctuations  measured  at  different  receivers 
for  the  same  event  is  another  characteristic  of  the  scattering  studied  here 
using  the  "edited"  data  set.  In  order  that  Chernov's  theory  be  applicable,  the 
separation  between  any  given  pair  of  instruments  is  projected  along  two  direc¬ 
tions:  the  "longitudinal"  direction  which  is  parallel  to  the  direction  of  pro¬ 
pagation  of  the  incident  wave  and  the  "transverse"  direction  which  runs  per¬ 
pendicular  to  the  longitudinal.  The  separation  between  any  two  given  receivers 
measured  with  this  coordinate  system  will  depend  upon  the  azimuth  assumed  for 
the  incident  wave  and,  therefore,  will  differ  among  events.  Now  consider  these 
separations  to  be  measured  on  a  grid  consisting  of  overlapping  squares  of  side 
20  km  in  the  longitudinal  and  transverse  directions,  i.e.,  a  grid  made  up  of 
blocks  defined  by 

x  <  longitudinal  separation  <  x 

li  -  2i  (21) 

xli  =  °»  5’  10«  90  ^  x2i  =  20,  25,  30,  ...,  110  km 

y1±  <  transverse  separation  (22) 

yli  =  °*  5»  10»  *•*»  90  ^  ?2i  =  20’  25’  30»  110  ^ 
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The  blocks  are  defined  by  (21)  and  (22)  to  overlap  in  order  to  smooth  the 
data.  For  the  nine  events  of  the  edited  data  set,  the  number  of  pairs  of 
instruments  with  separations  that  fall  into  each  of  the  blocks  defined  by  (21) 
and  (22)  are  shown  in  Figure  9.  If  (i,j)^  denotes  the  indices  of  a  pair  of 
instruments  which  for  the  kth  event  have  a  separation  which  lies  within  a 
given  block  and  if  denotes  the  set  of  all  such  pairs,  then  the  correlation 
coefficient  of  amplitude  fluctuation  for  the  block  is  given  by 


9 

*b  =  2 


k=l  S, 


B  B  /(  2  (  E  bJ 

J  k-1  all  i 


all  j 


B*)>  1/2 

J  « 


(23) 


A  similar  expression  holds  for  the  correlation  coefficient  of  phase  fluctua¬ 
tions  R  . 
s 

The  results  of  evaluating  R^  and  Rg,  using  the  edited  data  set  (with  its 
two  sets  of  phase  residuals),  are  shown  in  Figures  10,  11,  and  12.  The  extreme 
left-hand  column  and  the  uppermost  row  in  these  figures,  corresponding  to  pure 
transverse  and  pure  longitudinal  separation  of  the  instruments,  should  be 
compared  to  theoretical  results  derived  in  Rivers  and  von  Seggern  (1978).  An 
idea  of  the  reliability  of  the  correlation  found  for  each  block  is  found  in 
Figure  9,  which  shows  the  number  of  observations  used  for  each  measurement. 

Figure  10  shows  that  distinct  regions  of  positive  and  negative  amplitude 
correlation  exist.  The  nodal  curves  separating  these  regions  are  roughly 
linear  and  parallel  to  the  longitudinal  and  transverse  axes.  There  are  nodes 
at  transverse  separations  of  25  and  75  km  and  at  a  longitudinal  separation  of 
70  km.  The  amplitude  fluctuations  are  correlated  over  a  longer  distance 
(by  a  factor  of  almost  three)  in  the  longitudinal  direction  than  in  the  trans¬ 
verse  direction,  a  fact  at  least  in  qualitative  agreement  both  with  Chernov's 
theory  and  with  the  empirical  results  of  Mack’s  (1972)  study  of  Rayleigh-wave 
coherence. 

The  correlation  coefficient  for  phase  fluctuations  calculated  for  phase 
residuals  determined  on  the  basis  of  an  assumed  constant  C  for  all  events 
(Figure  11)  does  not  reflect  so  clear  a  spatial  pattern  as  does  the  amplitdue 
correlation  coefficient.  However,  the  pattern  becomes  clearer  when  the  resi¬ 
duals  used  are  those  which  were  calculated  using  a  different  value  of  C  for 
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TRANSVERSE  SEPARATION  (km) 


LONGITUDINAL  SEPARATION  (km) 

Figure  9.  Number  of  pairs  of  receivers  at  NORSAR  which  were  used  in  the 

edited  data  set  and  which  have  a  given  separation  (within  over¬ 
lapping  separation  intervals  of  size  20  km  X  20  km).  Note  that  the 
interval  size  is  actually  smaller  than  20  km  X  20  km  for  transverse 
and/or  longitudinal  separations  of  less  than  20  km. 
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Figure  10.  Correlation  (in  percent)  of  logarithmic  amplitude  fluctuations  at 
NORSAR  as  a  function  of  the  separation  between  the  receivers.  The 
amplitude  fluctuations  were  calculated  using  the  edited  data  set. 
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Figure  11.  Correlation  (in  percent)  of  spectral  phase  delays  at  NORSAR  as  a 
function  of  the  separation  between  the  receivers.  The  phase 
delays  calculated  using  the  edited  data  set  and  an  assumed  phase 
velocity  of  3.641  km/sec. 
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Figure  12.  Correlation  (in  percent)  of  spectral  phase  delays  at  NORSAR  as  a 
function  of  the  separation  between  the  receivers.  The  phase 
delays  were  calculated  using  the  edited  data  set  and  a  different 
value  of  the  phase  velocity  for  each  event. 
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each  event  (Figure  12).  For  this  second  set  of  phase  residuals,  R^  has  nodes 
at  transverse  separatic  ->  of  15  and  55  km  and  at  a  longitudinal  separation 
of  40  km.  Thus,  the  phase  fluctuations,  as  well  as  the  amplitude  fluctuations, 
are  correlated  over  nearly  three  times  as  great  a  distance  in  the  longitudinal 
direction  as  in  the  transverse  direction.  However,  the  nodes  for  R^  and  R 
occur  at  different  distances. 

Comparing  these  results  with  predictions  of  the  modified  Chernov  theory 
(Rivers  and  von  Seggern  (1978),  Tables  I-IV;  Figures  6-13)  reveals  certain 
important  results  and  discrepancies.  First,  the  theoretical  amplitude  corre¬ 
lation  in  the  transverse  direction  is  less  oscillatory  than  is  the  observed 
R^,  and  the  theoretical  value  does  not  oscillate  at  all  in  the  longitudinal 
direction.  Second,  the  theoretical  Rg  takes  on  no  negative  values  for  sepa¬ 
ration  in  either  direction  while  the  observations  do.  Reasons  for  these  dis¬ 
crepancies  may  lie,  at  least  in  part,  in  the  fact  that  the  values  of  R^  and  R 
for  large  separations  in  either  direction  are  not  particularly  reliable, 
because  they  are  based  on  few  observations.  Further  discrepancy  may  be  due 
to  the  constraint  that,  on  account  of  the  methods  used  in  calculating  them,  about 
half  of  the  residuals  must  be  positive  and  about  half  must  be  negative.  This 
result  again  reveals  the  problems  in  determining  the  correct  average  amplitude 
(or  phase)  from  which  the  fluctuations  are  to  be  calculated.  For  example,  if  an 
event's  average  amplitude  were  significantly  lower  in  the  vicinity  of  NORSAR 
than  was  the  average  of  the  values  measured  at  NORSAR  alone,  then  all  twenty-two 
amplitude  residuals  would  actually  be  positive,  and  R^  would  be  positive  for 
all  blocks  in  the  grid.  In  general,  then,  these  discrepancies  restrict  the 
comparison  with  the  theoretical  results  to  a  broad  discussion  of  the  amplitude 
correlation  coefficient  R^  in  the  transverse  direction. 

Figures  10  and  13  of  Rivers  and  von  Seggern  (1978)  reveal  that  for  large 

values  of  the  wave  parameter  D,  R^  has  a  node  in  the  transverse  direction  at 

£/ a  2  1,  where  £  is  the  transverse  separation  between  instruments  and  where 

a  is  the  correlation  distance.  Specifically,  Figure  13  of  that  report  shows 

a  node  at  l  =  1.15a  for  D  =  24.4  and  at  l  =  1.30a  for  D  =  85.6.  (The  difficulty 

in  calculating  D  from  our  data  sets  was  explained  above.)  Temporarily  consider 

the  value  D  =43,  determined  using  seven  events  from  edited  data  sets 

curve 
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that  were  analyzed  by  considering  three  distinct  periods,  a  common  azimuth  for 

the  waves  of  each  period,  and  a  different  phase  velocity  for  each  event  (Table 

XVIII).  For  D  ~  40,  the  node  occurs  at  about  SL  ~  1.2a,  so  a_  ~  (25km)/1.2  ~  21 

km.  For  waves  of  period  19.7  sec  and  phase  velocity  3.64  km/sec,  the  wave- 

number  k  =  0.0876,  so  ka  ~  1.8.  The  condition  that  ka  >>  1,  an  assumption 

2 

made  in  the  Chernov  theory,  is,  therefore,  not  valid.  Since  D  =  4L/ka  ,  the 

length  L  of  the  path  through  the  scattering  medium  is  found  to  be  roughly 

380  km,  a  value  that  is  much  less  than  the  actual  epicenter-to-NORSAR  distance 

for  any  of  the  events.  This  result  indicates  that  all  the  waves  incident  at 

NORSAR  have  been  scattered  once  within  380  km  of  the  array.  Since  details  of 

the  propagation  along  the  previous  2600  km  of  the  path  are  lost  by  this  scattering, 

we  cannct  say  whether  all  the  scattered  waves  were  part  of  the  same  wavefront 

at  a  distance  of  380  km,  as  in  the  P-wave  model  of  Aki  (1973)  and  Capon  (1974) 

or  whether  the  wavefield  at  NORSAR  is  the  final  result  of  the  multiple  scattering 

of  several  different  rays  which  diverge  over  a  longer  distance,  possibly  even 

back  to  the  source.  Had  the  opposite  approach  been  adopted  of  assuming  values 

of  about  21  km  for  the  correlation  distance  and  3000  km  for  the  path  length,  a 

value  of  about  310  would  have  been  obtained  for  the  wave  parameter.  Note  that 

the  many  disallowed  values  of  B^/S^  >  1  and  p  <  0  suggest  a  large  value  of  D 

since  many  of  these  disallowed  values  at  least  lie  within  the  confidence  limits 

surrounding  the  tail  end  of  the  bVs7  versus  p  curve,  for  which  D  -+  °°.  It  is 

unlikely,  however,  that  the  Born  approximation  could  hold  for  a  distance  L  of 

3000  km,  since  it  is  strictly  valid  only  for  single,  not  multiple,  scattering. 

Inability  to  solve  for  accurate  values  of  D  and  L  renders  it  impossible  to  solve 

accurately  for  the  remaining  parameter  in  the  Chernov  theory,  the  mean-square 

fluctuation  in  refractive  index  It  is  possible,  however,  to  calculate  an 

approximate  value  of  by  using  the  data  set  mentioned  above,  for  which  Dcurve 

=  43.  For  this  data  set  the  mean-square  amplitude  and  phase  fluctuation  are 
2  2 

B  =  0.0848  and  S  =  0.140.  The  value  of  may  be  calculated  by  using  either 

2  2 

of  these  values  and  the  formulas  for  B  or  S  given  in  Table  I.  It  is  assumed 
that  l)  =  43,  k  =  0.0876,  a  =  21^  km,  and  L  =  380  km.  Substitution  of  the  assumed 
values  into  the  formulas  for  B  yields  an  rms  fluctuation  of  /  ^'2'  =  4.46% 
and  substitution  into  the  formulas  for  yields  =  4.60%.  Although  these 

two  values  of  the  refractive  index  fluctuation  are  in  good  agreement,  their 
derivation  is  hardly  convincing,  so  they  should  be  regarded  only  as  rough  estimates. 
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RAYLEIGH-WAVE  SCATTERING  OVER  A  LINEAR  ARRAY  IN  THE  SOUTHWESTERN 

UNITED  STATES 


Data 

7 a  January  and  February  1962  eleven  LRSM  stations,  set  up  to  monitor 
GNOME,  an  explosion  in  southern  New  Mexico,  (Figure  1),  were  operative  in  the 
Southwestern  United  States.  These  eleven  stations  were  selected  to  form  an 
approximately  linear  array  extending  over  a  distance  of  almost  1100  km  from 
northern  Arizona  to  the  Big  Bend  region  of  Texas.  Analysis  of  the  amplitude 
and  phase  fluctuations  for  events  recorded  by  this  array  will  provide  insights 
into  the  effects  of  Rayleigh-wave  scattering  over  long  distances. 

The  two  events  chosen  for  this  study  were  an  earthquake  near  Nicaragua 
and  one  in  the  Kermadec  Islands  region  (Table  XIX).  These  events  were  selected 
so  that  incident  waves  would  be  both  longitudinal  and  transverse  to  the  separation 
between  the  stations  in  the  array.  It  would  of  course  have  been  helpful  to 
have  had  a  larger  data  base  than  two  events,  but  during  the  period  of  operation 
of  the  array  no  other  suitably  large  events  occurred  for  which  the  array-to- 
epicenter  azimuth  was  either  nearly  longitudinal  or  transverse.  Although  the 
ray  path  for  waves  from  the  Kermadec  event  crosses  an  oceanic-continental 
boundary  and  although  the  path  for  the  Nicaragua  event  runs  along  the  margin 
of  the  Pacific  Ocean  and  through  the  Sierra  Madre  mountain  range,  results  of 
a  ray-tracing  program  (Sobel  and  von  Seggern,  1978),  shown  in  Figure  13  and 
14,  reveal  no  significant  effects  at  the  array  due  to  systematic  multipathing 
(i.e.,  non-random  scattering). 

Computational  Methods 

Spectra  were  calculated  for  both  events  with  the  method  used  for  the 
NORSAR  data:  the  logarithmic  amplitudes  were  smoothed  over  the  band  of  periods 
18.3  to  21.3  sec,  and  the  phase  delays  measured  at  the  distinct  period  19.7 
sec.  However,  unlike  the  NORSAR  study,  where  the  same  time  window  was  anal¬ 
yzed  for  each  instrument,  the  large  dimensions  of  the  array  required  choosing 
a  time  window  for  analysis  that  was  different  for  each  station,  so  at  each 
station  the  maximum  amplitude  for  20-sec  Rayleigh  waves  would  arrive  centered 

Sobel,  P.  A.  and  D.  H.  von  Seggern  (1978).  Applications  of  surface-wave  ray 

tracing,  Bull.  Seism.  Soc.  Am.,  in  press. 
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Amplitudes  smoothed  over  0.0469  -  0.0547  Hz  band. 
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within  the  chosen  window.  After  visually  inspecting  the  seismograms,  we 
decided  to  select  a  time  window  that  propagated,  along  with  the  wave  train, 
across  the  array  at  a  velocity  of  3.87  km/sec.  If  the  true  phase  velocity 
had  been  equal  to  this  value,  the  phase  delays  calculated  by  the  spectral 
analysis  would  have  been  random  (as  it  was,  the  phase  delays  exhibited  a 
slight  linear  trend).  The  phase  fluctuations  were  assumed  equal  to  the  dif¬ 
ference  between  the  value  corresponding  to  this  linear  trend,  as  determined 
by  a  least-squares  fit,  and  the  observed  spectral  phases.  The  slope  of  the 
trend  removed  was  then  used  to  calculate  a  correction  to  the  assumed  value 
of  the  phase  velocity. 

Results 

The  amplitude  and  phase  fluctuations  for  the  two  events  are  shown  in 

Table  XX.  The  column  headed  "corrected"  contains  spectral  phases  with  multi- 

I les  of  360  degrees  added  or  subtracted,  as  necessary,  to  minimize  phase 

fluctuations  while  at  the  same  time  yielding  a  physically  realistic  value  of 

the  phase  velocity.  Note  that  other  choices  of  the  "corrected"  phase  would 

2  2 

have  resulted  in  different  values  of  C,  B  /S  ,  and  p.  While  for  the  Kermadec 
earthquake  certain  of  these  other  choices  cannot  be  completely  ruled  out,  for 
the  Nicaragua  earthquake  the  choice  shown  in  Table  XX  appears  valid. 

For  the  Kermadec  earthquake,  the  phase  velocity  C  =  3.37  km/sec,  the 

2  2 

mean-square  amplitude  and  phase  fluctuations  are  B  =  0.323  and  S  =  2.666, 

and  the  correlation  between  amplitude  and  phase  fluctuations  is  p  =  0.740. 

2  2 

For  the  Nicaragua  earthquake,  C  =  3,29  km/sec,  B  =  0.273,  S  =  0.410,  and 

~2 

P  =  0.056.  For  fluctuations  this  large,  particularly  the  case  of  S  for  the 

Kermadec  earthquake,  Chernov’ s  assumption  of  small  perturbations  is  clearly 

2  2 

violated.  The  ratio  B  /S  <  1  for  both  events,  but  for  the  Kermadec  earth¬ 
quake  p  >  0.5  and  for  the  Nicaragua  earthquake  p<0,  so  the  wave  parameter  D 

is  undefined  in  both  cases . _  However,  again  note  that  the  error  bars  are  quite 

2  2 

large  for  values  of  B  /S  and  p  determined  from  only  ten  observations.  (On 
the  day  of  each  event  only  ten  of  the  eleven  stations  were  operative.) 

For  the  spatial  autocorrelation  of  amplitude  fluctuation  R^,  note  that 
each  event  has  (^)  =  45  different  pairs  of  receivers  and  that  the  longitud¬ 
inal  or  transverse  separation  represented  by  each  pair  is  an  erratically 
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distributed  function  because  the  instruments  were  not  spaced  evenly  along  the 
array.  Therefore,  calculating  results  by  dividing  the  separations  into  seg¬ 
ments,  like  for  the  NORSAR  study,  will  fail  because  too  few  pairs  of  receivers 
exist  whose  separation  falls  within  a  given  segment  and  because  a  given  seg¬ 
ment  would  need  to  encompass  too  great  a  range  of  separations.  One  aji  hoc 
solution  to  this  dilemma  was  a  procedure  where  the  data  were  smoothed  before 
computing  the  correlation.  An  expanded  "data  set"  was  formed  by  interpolating 
the  values  of  the  logarithmic  amplitudes  at  10  km  increments  along  the  length 
of  the  array.  The  spatial  correlation  coefficient  was  then  calculated  for 
the  even  spacing  of  receiver  separations  this  expanded  data  set  represented. 

The  validity  of  this  interpolation  process  can  be  questioned,  however,  because 
in  practice  the  amplitudes  that  would  have  been  observed  at  the  evenly-spaced 
points  would  not  equal  the  interpolated  values.  Justification  for  this  proce¬ 
dure,  admittedly  questionable,  is  based  upon  the  NORSAR  study  which  showed 
that  the  rms  amplitude  fluctuation  over  the  separation  of  100  km  is  typically 
0.25  natural  log  unit  while  the  amplitude  fluctuations  for  the  linear  array 
are  more  than  twice  as  great.  Thus,  the  fluctuations  about  the  values  inter¬ 
polated  between  consecutive  stations  typically  100  km  apart  are  small  when 
compared  to  the  fluctuations  measured  between  stations  in  the  array  that  are 
several  hundred  km  apart.  Note  that  when  examining  the  results  of  this  pro¬ 
cedure  the  linear  interpolation  imposes  a  priori  the  constraint  of  strong 
correlation  over  small  separations.  Therefore,  values  found  for  separations 
of  less  than  about  100  km,  which  represents  the  true  resolution  of  the  linear 
array,  are  rather  tenuous. 

Figures  15  and  16  illustrate  the  amplitudes  measured  at  the  different 
stations  for  the  two  events.  The  value  of  R^,  determined  for  the  interpolated 
(dashed)  curves,  is  shown  in  Figures  17  and  18.  Again,  assuming  that  the  cor¬ 
relation  distance  ,a  of  the  scattering  medium  is  roughly  equal  to  the  distance 

where  R^C^-)  falls  to  zero,  this  value  is  estimated  to  be  about  125  km.  These 
figures  also  show  that  the  correlation  extends  in  the  longitudinal  direction 
about  1.8  times  farther  than  in  the  transverse  direction.  The  effect  of  atten¬ 
uation,  previously  ignored,  can  be  accounted  for  by  adding  a  correction  factor 
of  1.66  log  A  to  the  logarithmic  amplitudes,  where  A  is  the  epicenter-to— sta¬ 
tion  distance  in  degrees  (Figure  19).  This  term  is  nearly  the  same  for  all 
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Logarithmic  ampltudes  for  LR20,  not  corrected  for  attenuation  or 
station  effects,  for  the  Kermadec  earthquake. 


Logarithmic  amplitudes  for  LR20,  not  corrected  for  attenuation  or 
station  effects,  for  the  Nicaragua  earthquake. 


TION 


SEPARATION  (km) 

Figure  18.  Spatial  correlation  of  dashed  curve  in  Figure  16. 


6  80 


Logarithmic  amplitudes  for  LR20,  corrected  for  attentuation  but 
not  for  station  effects,  for  the  Nicaragua  earthquake. 


stations  for  the  Kermadec  earthquake,  so  only  its  effect  upon  the  Nicaragua 
earthquake  data  will  be  considered.  Figure  20  shows  that  this  correction 
extends  the  correlation  in  the  longitudinal  direction  to  separations  of 
about  265  km. 

Possible  amplitude  biases  at  the  eleven  stations  in  the  array  were  calcu 
lated  by  measuring  visual  magnitudes  at  each  of  the  stations  for  a  suite  of 
ten  events  (Table  XIX).  These  events  were  selected  for  their  azimuth  distri¬ 
bution,  so  the  effects  from  longitudinal  or  transverse  orientation  would  not 
be  subtracted  out  of  the  data  as  station  effects.  The  residuals  found  for 
each  station  from  these  ten  events  (Table  XXI)  were  removed  from  the  logarith 
mic  amplitudes  for  both  the  Kermadec  and  Nicaragua  events,  and  was  calcu¬ 
lated  again  (shown  in  Figures  21-24).  This  correction  effectively  decreased 
the  correlation  distance  for  the  Nicaragua  event  and  increased  it  for  the 
Kermadec  event.  When  both  the  attenuation  and  station-effect  corrections  are 
entered,  as  in  Figures  25  and  26,  the  correlation  extends  about  twice  as  far 
in  the  longitudinal  direction  as  in  the  transverse  direction. 

The  magnitude  of  the  phase  fluctuations  made  calculating  Rg  impossible. 
Since  many  of  the  spectral  phases  for  the  Kermadec  earthquake  could  not  be 
determined  uniquely  to  within  360  degrees,  the  reliability  of  the  phase  data 
is  highly  suspect  in  many  cases. 
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Figure  23.  Spatial  correlation  of  dashed  curve  in  Figure  21. 
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Figure  25.  Logarithmic  amplitudes  for  LR20,  corrected  for  attenuation  and 
station  effects,  for  the  Nicaragua  earthquake. 


IMPLICATIONS  OF  WAVE  SCATTERING  IN  THE  USE  OF  RAYLEIGH  WAVES 
TO  DETERMINE  SEISMIC  SOURCE  PARAMETERS  AND  EARTH  PROPERTIES 

Magnitude  Estimates  from  Long-Period  Rayleigh  Waves 

Event  Magnitude  (Mg)  from  Rayleigh  waves  is  commonlv  computed  on  the  basis 
of  amplitude  measurements  at  periods  near  20  sec.  Previously,  this  report  has 
dealt  with  amplitude  only  in  the  spectral  domain,  but  we  will  now  show  that  the 
spectral  amplitude  scatter  is  indicative  of  the  time-domain  scatter.  For  the 
events  studied  earlier  at  NORSAR  (Table  II),  maximum  seismogram  amplitudes  for 
periods  near  20  sec  were  measured  and  reduced  by  the  system  response  at  the 
exact  measured  period.  Because  the  propagation  distance  over  NORSAR  is  small, 
no  correction  is  required  for  differences  in  epicentral  distance  for  the  22 
individual  sensors,  and  the  variance  of  over  the  sensors  can  be  computed 
directly;  the  results  are  listed  in  Table  XXII  for  the  10  events.  The  average 
standard  deviation  of  M^  is  roughly  .09.  For  comparison,  the  results  using 
average  spectral  amplitudes  over  a  narrow  (.0488  -  .0508  Hz)  and  wide  (.0315  - 
.0908)  band  are  also  listed;  the  narrow  and  wide  bands  comprise  2  and  11 
points,  respectively,  in  the  257-point  spectra  (512  seconds  sampled  at  1 
s/sec).  Narrow  band  (.0468  -.0547  Hz)  estimates  were  used  in  previous  sections 
of  this  report.  Table  XXII  indicates  that  the  scatter  in  spectral  amplitude 
measurements  is  nearly  the  same  as  the  seismogram  measurements  and  also  indi¬ 
cates  that  a  wider  bandwidth  estimate  will  reduce  the  amplitude  variance  over 
NORSAR.  Using  the  results  given  in  Table  III  for  average  sensor  bias  at 
NORSAR,  a  standard  deviation  of  .04  Mg  units  results  from  only  the  purely 
local  effects  of  the  sensors,  neglecting  actual  wave  scattering.  The  sensor 
biases  are  based  on  only  10  events  and  would  probably  be  less  with  a  larger 
sample.  In  fact,  on  physical  grounds,  NORSAR  sensor  effects  are  expected  to 
be  negligible  for  20-sec  waves  where  the  wavelength  is  almost  as  long  as  the 
entire  array  aperture.  Therefore,  the  "mean"  o  values  in  Table  XXII  should 
be  reduced  only  slightly,  if  at  all,  to  give  the  true  Rayleigh-wave  scattering 
effects  over  NORSAR. 

Chang  and  von  Seggern  (1977)  studied  m^  scatter  across  LASA,  which  has 
13  subarrays  and  an  aperture  of  roughly  55  km.  They  reported  a  mean  o  for 
of  .15  based  on  a  large  event  sample.  Local  effects  at  the  subarrays,  they 
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stated,  accounted  for  only  a  small  fraction  of  this  fluctuation.  Because  the 
NORSAR  array  spans  over  100  km,  a  better  comparison  with  m^  scatter  can  be 
made  if  only  amplitudes  from  a  reduced  NORSAR  are  considered.  An  array  formed 
by  the  inner  seven  sensors  would  have  an  aperture  of  roughly  50  km,  nearly 
equivalent  to  LASA.  When  standard  deviations  of  spectral  amplitudes  over  these 
seven  sensors  were  computed  for  each  of  the  ten  events,  they  equalled  roughly 
one-half  the  values  for  the  full  NORSAR  array,  in  the  range  .04  -  .05.  When 
compared  to  .15  m^  scatter  over  LASA,  a  threefold  decrease  in  magnitude  scatter 
apparently  results  for  Mg  compared  to  that  for  m^  over  an  equivalent  array 
aperture.  Using  the  seismic  reciprocity  theorem,  Chang  and  von  Seggern 
argued  that  the  observed  level  of  receiver  fluctuation  for  might  also  apply 
to  expected  fluctuations  at  one  distant  receiver  from  equivalent  sources 
spread  over  an  area  with  dimensions  and  physical  properties  similar  to  LASA's. 
In  this  study  the  identical  argument  is  made  for  Mg  scatter,  with  the  infer¬ 
ence  that  Mg  is  a  more  stable  measure  of  source  strength  as  the  source  loca¬ 
tion  is  varied.  Therefore,  calibration  of  a  source  area  will  hold  to  larger 
inter-source  distances  for  Mg  than  for  m^.  Furthermore,  results  already  stated 
in  this  report  clearly  show  that  this  calibration  will  hold  to  significantly 
larger  inter-source  distances  measured  along  the  great-circle  path  to  the 
receiver  than  transverse  to  it. 

Seismograms  for  the  events  recorded  in  the  Southwestern  United  States, 
studied  earlier  in  this  report  in  the  spectral  domain,  provide  additonal 
insight  into  Mg  scatter.  Computations  identical  to  those  for  the  NORSAR 
events  in  Table  XXII  are  shown  for  these  two  events  at  the  bottom  of  the 
same  table.  The  result  revealed  that  1)  the  Mg  scatter  over  the  Southwestern 
U.  S.  array  is  roughly  twice  that  of  the  NORSAR  array;  2)  the  wide-band  spec¬ 
tral  measure  of  M  is  more  stable  than  either  analyst  M  or  a  narrow-band 
s  s 

spectral  estimate  near  f  =  .05  Hz;  3)  the  gross  signal  fluctuation,  measured 
in  either  the  time  or  frequency  domain,  is  nearly  the  same  for  sensors  aligned 
along  the  great-circle  path  as  for  sensors  perpendicular  to  it.  Also  in  this 
report  we  stated  that  for  the  Southwestern  U.  S.  array  the  signal  correlation 
was  definitely  higher  over  a  longer  space  dimension  for  the  Nicaragua  earth¬ 
quake  than  for  the  Kermadec  earthquake.  However,  because  the  dimension  of  the 
array  (~  1100  km)  is  large  compared  to  the  correlation  length,  the  overall 


level  of  the  observed  amplitude  fluctuations  does  not  differ  with  the  angle  of 
approach. 

Attenuation  Estimates  from  Long-Period  Rayleigh  Waves 

By  expressing  the  observed  Rayleigh-wave  spectrum  as: 

A(f )  =  S(f)  *G(A)  •e“a(f)r 


where  S  is  the  source  spectrum,  G  is  the  geometrical  spreading  factor,  and  r 
is  the  distance,  it  is  possible  to  take  the  ratio  of  the  spectra  at  two  sta¬ 
tions  on  a  great-circle  path  as 


Ax(f) 

A^TfT 


G(^) 


-Of  (f  ) 


(ri  ■  r2) 


By  absorbing  the  known  G  factors  into  the  amplitude,  we  find  for  the  attenuation 
coefficient 


ot(f)  =  [logA1(f)  -  logA2(f)]/(r2  -r2) 

For  20  sec  Rayleigh  waves,  a  typical  a  is  .0003  km  1  (Mitchell  et  al.,  1976). 
For  this  value  as  .0002  and  .0004,  Figure  27  shows  the  term  log  A^(f)  -  log 
A2(f)  versus  r2  -  r^.  Following  the  authors’  observations,  the  measured  dif¬ 
ference  of  log  A^(f)  -  log  A2(f)  can  have  several  tenths  error  due  to  random 
scattering  effects.  Therefore,  apparently  at  a  station  separation  of  even 
5000  km,  a  true  a  of  .0003  km  ^  could  be  estimated  at  anywhere  from  .0002  to 
.0004.  This  range  is  actually  wider  than  necessary  because  it  is  based  mei -  y 
on  scattering  data  for  one  frequency  while  attenuation  measurements  are 
usually  spread  over  a  fairly  wide  bandwidth  and  smoothed  in  some  manner. 

Phase-Velocity  Estimates  from  Long-Period  Rayleigh  Waves 

Data  in  Tables  X  and  XI  indicate  that  a  typical  standard  deviation  of  the 
measured  phase  over  the  22  NORSAR  sensors  is  .03  circles  (12°)  for  Rayleigh 
waves  at  .05  Hz.  For  the  two  station  method  of  determining  phase  velocity, 
the  equation 


C  =  f(r2  -  ri)/(4>2  -  4^  ±  n) 


Mitchell,  B.  J.,  L.  W.  B.  Leite,  Y.  K.  Yu,  and  R.  B.  Herrmann  (1976).  Attenua¬ 
tion  of  Love  and  Rayleigh  waves  across  the  Pacific  at  periods  between 
15  and  110  seconds,  Bull.  Seism.  Soc.  Am.,  66,  1189-1202. 
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is  used  where  C  is  phase  velocity  in  km/sec,  f  is  frequency,  <j>  is  the  measured 
phase  in  circles,  and  n  is  the  proper  integer  to  be  adjusted  to  make  C  phys¬ 
ically  acceptable.  For  small  errors  A<j>  in  -  <j>^  due  to  scattering 

C  +  AC  =  f  (r2  -r1)/(<i)2  -  ^  +  A<j>1  in) 

Using  C  =  4.0  and  f  =  .05  Hz  (oceanic  structure),  the  difference  AC  is  ±  .04 
km/sec  for  A<p  =  ±  .03  circles  at  r2  -  r^  “  200  km  and  decreases  as  r2  -  r^ 
increases,  assuming  that  A(j>  remains  bounded  for  receiver  separations  much 
greater  than  the  correlation  distance.  Thus,  to  yield  a  general  accuracy  of 
± 1%  to  phase-velocity  measurements  with  fair  confidence,  a  two-station  separ¬ 
ation  of  only  few  hundred  km  would  be  sufficient  if  random  scattering  effects, 
like  those  discussed  here,  were  the  only  factors  in  measurement  error. 

Source  Mechanism  from  Surface  Waves 

The  inversion  of  surface-wave  spectra  for  obtaining  the  source  mechanism  of 
of  earthquakes  has  been  widely  reported  as  a  successful  method.  Recently, 
estimates  of  the  number  of  spectra  necessary  for  resolving  the  source  mechanism 
have  decreased;  in  fact,  Turnbull  (1976)  suggested  that  as  few  as  two  stations 
would  be  adequate  if  placed  favorably  in  relation  to  the  source  radiation  pat¬ 
tern.  However,  several  propagation  effects  seem  to  argue  against  Turnbull's 
view,  especially  laterally  varying  attenuation  properties  and  large-scale 
refraction  and  multipathing.  In  this  report,  the  data  showed  that  only  small- 
scale  scattering  affects  the  signals  recorded  across  the  Southwestern  U.S. 
array,  but  that  the  effect  was  still  considerable  on  amplitudes.  As  a  guide 
to  the  spectral  variation  expected  for  recordings  of  an  earthquake  solely 
because  of  this  effect,  the  LR  spectra  for  the  Kermadec  and  Nicaragua  events 
are  shown  in  Figure  28  and  29  respectively.  Note  the  acute  fluctuations  for 
various  frequencies,  some  of  which  might  be  interpreted  as  spectral  "holes" 
related  to  source  depth.  Although  the  corresponding  Love-wave  specta  were  not 
processed  here,  the  authors  presume  they  display  similar  behavior.  In  view  of 
the  variations  seen  for  these  two  suites  of  recordings,  accepting  that  LQ  and 
LR  for  as  few  as  two  stations  can  reliably  constrain  the  earthquake  source 
mechanism  is  difficult. 

Turnbull,  L.S.  (1976).  Determination  of  seismic  source  parameters  using  far- 

field  surface  wave  spectra,  Ph.D.  Thesis,  Pennsylvania  State  University 

University  Park,  Pennsylvania. 
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Figure  28.  Spectra  of  the  Nicaragua  earthquake  LR  recordings  at  a  linear 
array  of  LRSM  sites  in  the  Southwestern  United  States. 
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Figure  29.  Spectra  of  the  Kermadec  earthquake  LR  recordings  at  a  linear 
array  of  LRSM  sites  in  the  Southwestern  United  States. 
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CONCLUSIONS 


Both  amplitude  and  phase  fluctuations  of  20-sec  Rayleigh  waves  were  found 
to  correlate  over  a  longer  distance  in  the  direction  of  wave  propagation  than 
perpendicular  to  it.  In  addition,  the  report  revealed  that  amplitude  and  phase 
fluctuations  were  correlated  positively,  in  most  cases.  While  this  fact  might 
support  the  hypothesis  of  random  scattering  of  the  Rayleigh  waves  according 
to  a  modified  Chernov  theory,  the  hypothesis  of  multipath  arrivals  cannot  be 
rejected.  Data  available  for  this  study  was  insufficient  to  infer  the  exact 
nature  of  the  wavefield.  In  the  case  of  the  NORSAR  investigation,  the  array 
aperture  of  100  km  was  too  small  to  define  the  correlation  functions  accurately 
to  a  large  enough  spatial  dimension.  In  the  case  of  the  Southwestern  United 
States  linear  array,  the  spacing  (about  100  km)  was  too  large  to  construct 
confidently  the  spatial  correlation  function  and  the  entire  aperture  of  1000 
km  was  so  extended  that  large-scale  multipathing,  although  not  predicted  by 
the  ray-tracing  program,  might  have  affected  the  observations.  From  studies 
of  recordings  at  the  Southwestern  U.S.  array,  indications  are  that  the  cor¬ 
relation  length  of  the  inhomogeneities  causing  the  random  scattering,  assuming 
it  exists,  is  the  order  of  100  km  or  more.  The  smaller-scale  scattering 
observed  at  NORSAR  is  characterized  by  a  correlation  length  of  about  25  km. 

The  standard  deviation  of  .05  Hz  spectral  amplitudes  from  Rayleigh-waves 

recorded  on  the  NORSAR  array  of  instruments,  roughly  100  km  in  aperture,  is 
approximately  0.1  log  unit.  For  the  same  events  used  in  the  spectral  analysis, 
the  analyst-measured  Mg  values  also  had  a  standard  deviation  of  0.1  log  unit. 

This  figure  is  only  somewhat  less  than  the  0.15  standard  deviation  reported  by 
Chang  and  von  Seggern  (1977)  for  m^  values  across  the  50-km  aperture  of  LASA. 
However,  when  only  the  inner  seven  NORSAR  sensors  are  considered  (aperture  about 
50  km),  the  Mg  standard  deviation  drops  to  roughly  0.05.  Thus,  Mg  scatter  is 
significantly  less  than  m^  scatter  over  equivalent  array  dimensions,  and  by 
reciprocity,  Mg  would  provide  more  accurate  (1/3  the  standard  deviation)  rela¬ 
tive  yield  estimations  for  shots  within  a  test  site  than  would  m^. 

On  the  basis  of  the  observed  level  of  amplitude  fluctuations  reported  here, 
we  suggest  that  attenuation  measurements  for  Rayleigh  waves  can  give  highly 
variable  results  for  differing  sample  points  of  the  wavefield  and  that  the 
attenuation  coefficient  can  be  resolved  confidently  to  within  a  factor  of  two  only 


for  long  paths,  at  least  of  few  thousand  kilometers.  On  the  other  hand,  the 
associated  phase  fluctuations  are  relatively  insignificant  in  phase-velocity 
measurements;  for  20-sec  LR,  the  phase  velocity  can  be  resolved  to  within 
1%  over  short  (>  200-300  km)  paths  if  instruments  are  well  calibrated  and  no 
large-scale  multipathing  occurs. 

Spectra  of  LR  waves  are  considerably  affected  by  small-scale  scattering, 
and  their  variability  would  seem  to  preclude  accurate  source  mechanism  resolu¬ 
tion  from  a  small  number  of  seismograms. 
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APPENDIX  I 

THEORY  OF  WAVE  PROPAGATION  IN  A  TWO-DEMENS IONAL  RANDOM  MEDIUM 

In  a  companion  report  (Rivers  and  von  Seggern,  1979)  the  authors  investi¬ 
gated  fluctuations  in  the  amplitudes  and  phases  of  20-sec  Rayleigh  waves  caused 
by  random  irregularities  in  the  structure  of  the  crust  along  the  great-circle 
path  of  wave  propagation.  By  analyzing  these  fluctuations,  which  were  measured 
at  NORSAR  and  in  the  southwestern  United  States,  certain  parameters  can  be 
evaluated  that  characterize  the  degree  of  heterogeneity  in  the  structure  of 
that  segment  of  the  crust  through  which  the  waves  travelled.  In  order  to  perform 
this  analysis  it  is  necessary  to  have  a  theoretical  prediction  of  the  scattering 
phenomena  that  are  expected  when  an  elastic  wave  propagates  through  an  inhomo¬ 
geneous  medium  characterized  by  continuous,  but  random,  fluctuations  in  the  index 
of  refraction.  Chernov  (1962)  developed  such  a  theory  of  wave  propagation  in 
three  dimensions.  In  order  to  study  the  Rayleigh-wave  problem,  the  authors 
have  modified  Chernov's  theory  to  make  it  applicable  to  the  two-dimensional 
case.  This  memorandum  consists  of  the  derivation  of  certain  theoretical  results 
summarized  in  Table  I  of  Rivers  and  von  Seggern  (1979),  against  which  the  ob¬ 
served  Rayleigh-wave  amplitude  and  phase  fluctuations  are  to  be  compared. 

Although  the  discussion  presented  here  of  the  two-dimensional  random 
scattering  problem  is  intended  to  be  more  or  less  self-contained,  it  neces¬ 
sarily  requires  knowledge  of  Chernov's  (1962)  text,  which  it  follows  closely. 
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Amplitude  and  Phase  Fluctuations 


We  assume  that  the  random  medium  is  characterized  by  only  small  deviations 
of  the  density  and  acoustic  velocity  from  their  mean  values,  i.e. 

=  >  Ac  -  c-c.  «-c8 ,  (1) 

We  assume  further  that  a  plane  wave  having  acoustic  pressure  pg  and  travelling 
in  the  +x  direction, 

~  A»  J-'Mj.  -  vui]  (2) 

is  incident  upon  the  medium.  Each  element  of  the  medium  then  becomes  a  source 
of  secondary  scattered  waves  p^  which  obey  the  equation  (Chernov,  p.  43,  (26))* 

^  "  V  ‘  Co  4  A  J  A0»^WUt-Kx)J.  (3) 

Taking  the  coordinates  of  the  scattering  element  to  be  (£,n),  those  of  the 
observations  point  to  be  (x,  y),  and  their  separation  to  be  r, 

r  =  J(x-£)1  (*) 

we  may  solve  the  inhomogeneous  equation  (3)  by  the  method  of  Green's  functions, 
obtaining 


\  x  ac  ±JL  1!A£2  1  '‘■*1  .  <»,  v  -i* 
))  U*-  +  p.  Je 


J-n 


where  the  integration  extends  over  the  entire  region  from  which  waves  are  being 
scattered  to  the  observation  point.  The  index  of  refraction  of  the  medium  is 
defined  to  be 


-  (C/M 


so  by  substituting  (6)  into  (5),  dropping  the  time  dependence  e  ,  and  neg¬ 
lecting  the  scattering  caused  by  density  fluctuations  in  comparison  with  that 
caused  by  velocity  fluctuations,  we  obtain 


x  a,  y 

P'  s  X 


lie  *  Hi0  Ik*-)  pU,^)  ^  ^ 


For  waves  of  frequency  0.05  Hz  and  phase  velocity  3.6  km/sec,  the  wave  number 
K  is  0.087  km  ,  -and  the  asymptotic  expansion  of  the  Hankel  function, 


i  (  k  r  -  m  /m) 


H „  (kr)  ~  4  vxr  C 


*  All  references  are  to  the  text  of  Chernov  (1960). 
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is  valid  to  within  a  few  percent  for  scatterer-to-receiver  separations  r 
of  more  than  about  60  km.  The  equation  obeyed  by  the  scattered  waves  thus 
becomes 

r77  ((  ei*l»**i>**M 

P>  ))  47  p  u^)  i  r,  i-n 

This  equation  is  the  two-dimensional  analog  of  Chernov,  p.  59,  (67). 

If  the  sum  of  the  incident  and  scattered  waves  is  expressed  as 


p  =  A  e 


-v 


(10) 


and  if  the  amplitude  A  =  A^  +  AA  and  phase  <J>  =  q>  +  Aq>  of  the  resulting 
wave  differ  only  slightly  from  those  of  the  incident  wave,  i.e., 


A  A. 

A.  «  i  1  A  «  i-j 


it  follows  (cf. Chernov,  pp.  58-60)  that 


(11) 


AA.  _  \\  ft  vj  (kIt-(x*  c,)^  -t-  ^  )  p  (i,m)  d g,  A-n 

A  o  ' 


(12) 


and  that 


Af  =  JJj  Si  4=  0*4Klr-U-O}+  h)|<MkK 


(13) 


Equation  (11)  is  known  as  the  Born  approximation  in  scattering  theory. 
Chernov  uses  the  method  of  Rytov  to  show  that  (12)  and  (13)  hold  under  conditions 
less  rigid  than  those  imposed  by  (11),  namely  that  the  amplitude  fluctuation 
AA/A0  and  phase  fluctuation  A<£  be  small  over  the  distance  of  a  single  wave¬ 
length,  rather  than  over  the  entire  path  length  of  the  wave  through  the 
medium.  This  is  an  important  loosening  of  restrictions,  since  the  observed 
fluctuations  are  not  always  small.  In  the  event  that  the  amplitude  fluc¬ 
tuations  are  not  to  be  regarded  as  small  perturbations,  Rytov' s  method 
dictates  that  the  left  side  of  equation  (12)  be  replaced  by  the  logarithmic 
amplitude  fluctuation  £n  AA/Aq.  We  must  note,  however,  that  Rytov' s  method 
has  been  criticized  by  Aki  (1973),  so  we  should  be  cautious  about  the  appli¬ 
cability  of  the  theory  to  those  cases  for  which  (11)  is  not  satisfied. 

Introducing  the  dimensionless  variables 


(14) 


and  denoting  the  phase  fluctuation  by  S  and  the  logarithmic  amplitude  fluctua- 
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tion  by  B,  equations  (13)  and  (12)  become 


6  U',/) :  Tjlf  ^  77  t/,A (. r ' ■  U'-%')  +  "T )  p(  ^V)  Jt,'  jy. 


(15) 


(16) 


We  will  treat  only  the  case  in  which  the  medium  is  composed  of  large-scale 
inhomogeneities,  i.e.. 


Ka  •»  1 


(17) 


where  a  is  the  characteristic  dimension  of  the  scatterers.  For  this  case 
we  are  able  to  neglect  back-scattering  and  assume  that  all  secondary  waves 
Pl  are  emitted  into  a  cone  (actually  a  two-dimensional  sector)  which  has  an 
aperture  angle  of  the  order  of  1/Ka  and  which  is  centered  about  the  direction 
of  the  incident  wave  pg  (i.e.,  the  +x-axis) (Chernov,  pp.  56-57).  It  follows 
that  the  receiver  at  (x‘,y')  will  detect  ratiation  scattered  only  from  those 
inhomogeneities  located  inside  a  similar  cone  with  apex  at  the  receiver  and 
oriented  in  the  -x-direction.  Within  this  cone  we  may  approximate  (4)  by 


r's  U'-  £')+  t 


(18) 


and  substitute  (18)  into  the  arguments  of  the  trigonometric  functions  in  (15) 
and  (16).  Since  l/Zr5"  varies  much  more  slowly  than  do  the  trigonometric  func¬ 
tions,  we  may  make  the  rougher  approximation 


4?  ~ 


(19) 


in  these  same  equations.  Since  we  need  consider  only  scatterers  lying  between 
the  planes  V =  0  and  £'  =  x' ,  (15)  and  (16)  become 


C  io*  ^  +  h)  (20) 


We  emphasize  that  our  neglect  of  waves  not  scattered  forward  into  a  small  cone 
(the  Fresnel  approximation  of  diffraction  theory)  is  valid  only  if  (17)  applies. 
If,  as  previously  supposed,  K  =  0.09  km-l,  our  theory  will  not  be  valid  for 
media  composed  of  scatterers  for  which  the  characteristic  dimension  is  much  less 
than  100  km.  Since  we  wish  to  consider  the  effect  of  inhomogeneities  smaller 
than  this,  it  will  be  necessary  to  make  approximations,  wherever  possible,  which 
are  not  so  strongly  dependent  upon  the  condition  imposed  by  (17). 
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We  assume  that  the  receiver  is  located  at  (L,0)  where  L  is  the  distance 
through  the  medium  travelled  by  incident  wave.  We  now  write  (20)  and  (21)  as 


'K 


where  we  have  adopted  the  abbreviations 

A  *  .  t  fi  2.  H  \ 

A,  +  hJ 


(22) 

(23) 


(24) 


( 


.*  ol 


(£♦*) 


(25) 


and  where  L'  =  K  L. 

For  the  moment,  let  us  make  explicit  the  fact  that  we  have  oriented  the 
+x-axis  of  our  coordinate  system  along  the  source-to-receiver  line  by  expressing 
coordinates  in  this  system  as  (x,y,0),  where  Q  is  the  geometrical  azimuth  of 
the  source  as  measured  in  a  geographical  (i.e.,  fixed)  coordinate  system 
centered  at  the  receiver.  Considering  coordinate  systems  oriented  toward 
two  different  sources,  we  see  that  the  points  (x^.y^.Q-^)  and  (xj_,  y^.Q?)  are 
physically  distinct,  but  in  the  two  different  coordinate  systems  they  Bear 
an  identical  geometrical  relationship  to  the  source  and  the  receiver  (see 
Figure  1).  (We  assume  all  sources  to  be  at  the  same  distance,  L,  from  the 
receiver.)  We  now  introduce  the  spatial  correlation  function  of  the  index 
of  refraction 


Nns  fl'4,8)  (26) 

where  the  overbar  denotes  the  average  taken  over  all  receiver-to-source  azimuths 
6,  i.e.. 


/  an 

Nn.  =  Tt  \, 


(27) 


Due  to  the  random  arrangement  of  scatterers,  the  correlation  function  N12  depends 
only  on  the  separation  between  the  points  (x]_,yj_,0)  and(x2,y2>©)  and  not  on  the 
points  themselves,  so  we  may  write 


where  is  the  mean  square  fluctuation  of  the  refractive  index  (i.e.,  N^) 
and  where  N  is  the  correlation  coefficient.  Now  the  geometrical  functions 
A^  and  A2  in  (22)  and  (23)  depend  only  on  the  coordinates  of  the  points 
(L,0)  and  (£,n)  and  not  on  any  physical  property  corresponding  to  the  actual 
points  in  space  themselves  (such  as  jj),  hence  they  are  the  same  whether  eval¬ 
uated  at  the  point  (£,n,Qi)  or  at  the  point  (£,n»02)*  It  then  follows  that  by 
squaring  (22)  and  (23)  and  by  averaging  over  all  receiver-to-source  azimuths 
9,  we  find  the  mean  square  phase  and  amplitude  fluctuations  to  be  » 


(29) 


(30) 


where 


r'=  )\ 


(31) 


We  rewrite  (29)  and  (30)  in  terms  of  new  coordinates  n  and  y,  defined  by 


(32) 


)  =  '+o, 


(33) 


Our  expressions  then  become 


sl  =  L  L  D  A .  ( V  -  j'j- T)Nir')  iMy  A?.'  (34) 

s’*  pi-  Alu'-%;Jy+T)AaU'-(.')'f  h:  h:.  os) 

It  is  shown  in  Appendix  II  that  performing  the  integration  with  respect  to  y 
yields 

s^P  lA  L*L  1 1  a, U'.'-t.'-t; ,•*>)+  u;  (36) 

B1  '-f  ^  Si  A  iVA,Uv.'-^-6.'^)  +  A.,V\V-^W]bUr (37) 

which  we  write  as 
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S'*  *  i  p  d.-li) 

where  we  have  defined 


(38) 


(39) 


i,»m' 

i,= ( 


O  f  o 


A.lU,'-  S*l/n)  Nlr')^u;  h: 


Introducing  new  variables  x  and  5,  defined  by 


(AO) 


(Al) 


X  - 


(A3) 


(AA) 

and  allowing  the  integration  over  i  to  extend  from  -<»>  to  +«>  (which  we  may  do 
if  L  >>  a,  since  N(r)  =  0  for  r  >>  a),  we  have 


I,'  L  S.„  i.M  A, l  Ulj'h)  N(f')  JK  k%  An  (45) 

1^  S.M-*  A,Ul/-lx/n)Nlc')JhUJx  (A6) 

where  we  now  write 

r'- 

(A7) 


Since  the  integrand  of  (A5)  is  independent  of  x,  we  find 


Ii-L  N(r')  «H  <*£• 


(A8) 


Recalling  our  definition  of  A^»  (2A),  and  noting  that  r'  is  an  even  function 
of  both  £  and  n,  we  write  this  as 
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Hi/  (“  (°°  j_ 

)o  ).  jf 


aw  (u  +  n  )  Nlr')  i"h  iz. 


\a  *),  -T&  \./^  )  +  c»a  1  at  )1  Nlr'jA’M  ij. 


(49) 


We  now  assume  (Chernov,  pp.  7-11)  that  the  correlation  coefficient  of  the 
refractive  index  has  the  form  of  a  Gaussian  density,  i.e.. 


-f'/a'  -r'Vo'*- 

e  5  e 


(50) 


where,  as  before,  a  is  the  characteristic  dimension  of  the  scatterers  and 
where  a1  -  Ka.  It  follows  from  (47),  (49),  and  (50)  that  by  setting 

n2 

■XT’  “  q  we  obtain 


<U 


-  L 


•it 


(51) 


21 

where  C  =  a77.  We  may  evaluate  the  integral  in  (51)  approximately  by  noting 
although  therange  of  integration  extendsfrom£  =  0  to  £  =  +°°,  the  exp(-£2/a'2) 
term  causes  the  integrand  to  vanish  for  £  >>  a’,  so  the  effective  range  of  £ 
is  actually  of  the  order  of  a'.  Then  throughout  this  effective  range  of  £  we 
have  £  -  2£  2 

— j-2  S  —i  «  1,  where  we  have  used  the  restriction  (17)  that  a'  = 

Ka  >>  1.  We  thus  set  £  =  0  in  (51),  obtaining 


l^U-'Ce' 


-  4V  K  La.  (52) 

which  is  identical  to  Chernov's  for  the  analogous  three-dimensional  case 
(Chernov,  p.  73  and  p.  83). 

We  may  obtain  a  more  cumbersome,  but  exact,  representation  of  1^  by 
writing  (51)  as 


I, 


(53) 
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where  a  *  J 'z;2  +  1)^^  +  *;  and  3  =  +  1)^^  -  C*  This  reduces  to 


i,  -  V  So  c  "  l  Jc’  + 1  +  ^ 1 


-  V  i  '  T  i  _l_  A, 

<?  4  »  a  * 


where  Z  =  /cz  +  1.  Equation  (54)  has  the  form  of  a  tabulated  Laplace 

transform,  and  it  reduces  to 

I.  -  T  K„Jr  , 


L  o“ 

xT  e 


k 


We  note  that  the  asymptotic  expansion  of  the  modified  Bessel  function  K^Cx), 


J  it  -  x 

K.  (x)  ~  j  x  £ 


/  \ 
(  U  **  t  .•  ■) 


implies  that  (55)  reduces  to  (52)  if  Ka  >>  1.  We  may  thus  use  (52)  when 
(17)  is  valid,  but  we  should  use  (55)  when  considering  scatterers  of  size 
a  <  100  km. 

Substituting  from  equations  (47)  and  (50)  into  (46)  and  performing 
the  integration  over  £,  we  find 


li  =  ^  o'  ).»  A,  (3L'  -ix.-h)  e  °'1 


Reversing  the  o.'der  of  integration  and  noting  that  the  integrand  is  an  even 
function  of  n,  wt  have 


1.  -  )o' 


'  (  "  ( L 


—  •  <—+■  + -) 
•  -  ~  C1-'  V  >u  '  -  “  X  +  “>  / 


|V.  S  X 


e  Ji  x  J  -p 


Substituting  t  *  /4L’  -  4x,  this  becomes  , 

T ,  -  a  I  r  i  ^  TT  ev.  (t  +  0  4  +  J  ^  I'M 

-  WL'^l  C  "r  '-•  0‘  +  e'  JC 

--  H„'L'  5.”  1  L*  C  e  '~"~r  Jr]  7  U1  +  T)jt 

' —  So  { c~  -  l)  ~L  av.  U  4 

-  _  SLl  r  “  I  -w1  ,i  i  ( *  +  ) 

5  ).  U  -Or  c  n 


where  we  have  adopted  Chernov's  (p.  74)  definition  of  the  so-called  "wave 
parameter"  D, 


D  =  ^ 


HL 

Kov 


(60) 


Physically,  the  wave  parameter  represents  approximately  the  ratio  of  the 
size  of  the  first  Fresnel  zone  ttXL  to  the  scale  of  the  inhomogeneities. 

Making  the  change  of  variable  p  =  e  i7T/^  t  in  equation  (59),  we  find 

-  ,  lOy'  S  J_ 


1  %  z  5l 

_  _al2  f  •" 

5.  ' .  a 

o.'4" 

"  H  K  t 


W' 


A  <  \)  a 


\ 


-i)  17  e 


x  -I./1 

e  jp 


V 


'  ’  (61) 
This  last  integral  has  the  form  of  a  Laplace  transform,  and  it  reduces  to 


—  1 
4  ! 


JrvT  K1  La 
0 


(62) 


We  finally  obtain  the  desired  expressions  for  the  phase  and  amplitude 
fluctuations  by  substituting  the  equation  for  Ij_,  either  (52)  or  (55),  and 
those  for  I2,  (60)  and  (62),  into  equations  (38)  and  (39).  We  see  that  the 
form  of  these  expressions  depends  in  an  essential  way  on  the  value  of  the 
wave  parameter  D: 


for  D  >>  1:  I2  =  0 

_  —  J¥  —  , 

=  p’  -  T  ^  *  L  q 

for  D  <<  1:  I2  =  /r  K2La 

=  JtT  pK'Ujf7  =0. 


(63) 


These  limiting  values  are  the  same  as  those  for  the  three-dimensional  case 
treated  by  Chernov  (p.  83,  (159)  and  (160)).  We  note  that  use  of  the  exact 
value  of  1^,  given  by  (55),  would  have  yielded  a  negative  value  for  B2  in  the 
case  D  <<  1.  We  also  note,  however,  that  the  limit  K  -*•  00  (which  would  imply 
D  -*■  0)  must  not  be  approached  too  closely,  since  our  scattering  formulas 
are  invalid  in  the  very  high  frequency  case  corresponding  to  geometrical 
acoustics  (Chernov,  p.  56). 
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Correlation  of  the  Amplitudes  and  Phase  Fluctuations  at  the  Receiver 

Having  determined  the  magnitude  of  both  the  amplitude  and  phase 
fluctuations,  we  now  examine  the  correlation  between  them.  We  deter¬ 
mine  the  average  (over  receiver-to-event  azimuths)  of  the  product  of 
the  two  fluctuations  by  the  same  method  as  was  used  to  determine  the 
average  of  their  squares,  (34)  and  (35),  thereby  obtaining 


A, [  j +  T  JMU ')AnJty  <K'  l%[.{ 


It  is  shown  in  the  Appendix  II  that  performing  the  integration  with  respect 
to  y  yields 


Changing  variables  as  indicated  in  (43)  and  (44)  and,  as  before,  allowing 
the  integration  with  respect  to  £  to  extend  from  £  =  to  £  =  +»,  we  find 

B5  -  >  p  L  ,U  l  «k  lx 


A1(lL/-lx)',.)NU')l'n  lx. 


We  now  substitute  (47)  and  (50)  into  (66)  and  perform  the  integration 
over  £,  obtaining 


I'u.'-nx) 


’l  .  (  00  ( l 


4  a'  a 

*-  '  v  l  o 


5-  ~  lo  )hV.  -4X 

By  inspection  of  equations  (58)-(61),  we  see  that 


V  .  V 

I  -n  '  Tx 

I  ^ _  3l\  -  2£ 

4  Hl'-H*  +  1  j  e  °  Hx 


--  -  f 

ol3 

H 

K 

(“  /  -*-PP'L  A  J-  XlT/4 

»o  U  - 1 )  e  e  dp 

.  i 

0 

f  °° , 

t 

t  \  e  p- 

4 

L  1 

Op  -  1  )  p'  dp 

V 

n'» 

e 

=  p 

*T 

S.  ( 

Or..  0  (i  "  ^  )  Vx  d  0  , 

l  1  1 
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Equation  (68)  is  tabulated  Laplace  transform,  so  finally  we  have 

—  -Hi  7' a' 5  r  i  /  ! -  \  *  i 


fcs  -  - 


p-  Li-feUiWl 


♦>)*]. 


The  correlation  coefficient  for  amplitude  fluctuations  and  phase 
fluctuations  is  given  by 


/>-  Jfi’Ji'  .  (70) 

If  we  approximate  in  (38)  and  (39)  by  (52),  we  have 

s1  =  *T  ^  K^aL  I  I  +  "p"  (.  lu’  +  l  -  0  ']  (71) 

p  KlaL  ll  -  %  (72) 

and  (70)  reduces  to 

-1  ^[l-fpTTi 

r-  usi-ijkti1  .  <73) 

We  see  that  for  large  values  of  the  wave  parameter  (D  >>  1)  we  have  p  =  0,  and 
for  small  values  of  the  wave  parameter  (D  <<  1)  we  have  p  =  1/2.  For  the  three- 
dimensional  case  Chernov  (p.  89)  finds  corresponding  values 
of  0  and  /3/2  /2,  respectively. 


Longitudinal  Autocorrelation  of  Amplitude  and  Phase  Fluctuations 

We  now  examine  the  correlation  of  both  amplitude  and  phase  fluctuations 
which  are  measured  at  two  separate  receivers.  By  "longitudinal"  autocorrela¬ 
tion  we  mean  the  special  case  in  which  the  separation  of  the  two  receivers  is 
parallel  to  the  direction  of  propagation  of  the  incident  wave  (i.e.,  normal  to 
the  wavefront).  We  thus  take  the  coordinates  of  the  receivers  to  be  (L^.O)  and 
(L2,0).  Taking  the  product  of  the  fluctuations  at  the  two  receivers  and  aver¬ 
aging  over  azimuths  (changing  the  actual  location  in  space  of  the  geometrical 
point  (L2,0,9)  so  that  the  separation  is  longitudinal  for  each  azimuth),  we  find 


where  L]/  “  KL^  and  L2 '  *  KL2 ■ 


We  assume  that  the  effect  of  waves  scattered  in  the  region  between  the  two 
receivers  is  negligible,  since  such  waves  are  incident  upon  only  the  second 
receiver.  We  thus  set  the  upper  limit  on  the  integrals  over  £2'  to  be  Lj_'. 
Inspection  of  equations  (29)- (41)  enables  us  to  write 


Vi,  --  if 

(76) 

Fj  1  c j  -  1 

(77) 

where  we  have  now  defined 


t  ^  / 

1,5  L'C  aaiu nu') dr  <u,'  <k' 
i, --  C: S'? LI a, { v :  ♦  l; - ( £/ ♦ 1; ),-*,)  nj uo a t; <h' 


(78) 

(79) 


Making  the  substitutions  (43)  and  (44)  and  setting  the  limits  on  the  £ 
integral  equal  to  and  -H»,  we  find 

Vi.TX  Uu;-L;-si 


(80) 


where  AL  =  L2  -  Lj_.  Proceeding  as  in  equations  (49)- (53),  we  have 

**  ]=f  \jc\T  +  l]  *  U 

where  now  l  -  ^  1AL  #  ye  divide  (81)  into  two  integrals 

aT  *- 

I.*  IC V5*  *'!'«  +i “c  e'^  .0S«}(»2) 


(81) 


where  * 

I,*  7*  i  S 


2 (AL1  -  £)  2 (£  -  AL') 

a ' 2  and  £2  *  a' 2  .  We  now  write 

At')"  « 

JT? 

'v£ 


*0VX,)+ 


^  s:  u-  ^ -A°' + e-  ■■■ ta^  * »]  ix .  a4 « 


Lla'1  _  AX  («  _  3_il  t  _  , 

IT  <?  )<.«■’  4  (al'  e)  JtsT  U^+>  +  0  "  <K . 
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(83) 


Inspection  of  (83)  reveals  that  when  the  separation  between  receivers,  AL, 
is  small  compared  to  the  size  of  the  scatterers,  a,  we  may  neglect  the 
dependence  of  I^upon  AL,  whereupon  1^  takes  on  the  value  found  in  the 
single-receiver  case. 


I.S-TnK'La,  (84> 

We  may  treat  the  general  case  by  evaluating  the  integral  in  (83)  numerically. 

We  substitute  (43),  (44),  (47),  and  (50)  into  equation  (79)  and  once 
again  set  the  limits  on  the  C  integral  equal  to  -<*>  and  +»,  thereby  obtaining 


-  X  "  ^  J  (-00^-00  j  XT  ll_'  +  u'.'XX  I  UiJr  ill  l- 


(  -XX)  +  H 


At 


=  J7(l7;I77~  (Mi'll.; 


r)-sx  +  tje  *  JxJ-n. 


By  comparison  with  equations  (58)- (62)  we  see  that 

,  t  l  f  T|  ^  '■  1  j_  ■  I  i  m  \  q  -  “~r% 

{.  L  j  v.i  j  t1  1  +  1  )  At  J  y\  e 


J^T  a'5 


>m  -i]1  -UoiTT  -0 1 } 


where  now 


t),  '  a'x  (  L  (  +  L  j ) 

k*  -7-  At'. 


For  the  case  in  which  the  receiver  separation,  AL,  is  small  compared  with 
the  source-to-receiver  distance,  L  s  L^  .  Ij,  the  expression  for  l£  reduces 
to  the  constant  value  found  for  the  single-receiver  case, 


-J  ITT  yOj-a 


K  \  \ —  11 

o  LJl/%1  -'J 


The  correlation  coefficient  for  phase  fluctuations  as  a  function  of 
receiver  separation  in  the  longitudinal  direction  is  given  by 


and  that  for  amplitude  fluctuations  is 

Rb,  long  “  p 

1,(01)-  lx[b'-) 

l.W'IJo)  (91) 


Transverse  Autocorrelation  of  Amplitude  and  Phase  Fluctuations 


We  now  examine  the  case  in  which  two  receivers  are  separated  in  the 
direction  "transverse"  to  the  direction  of  the  incident  wave  (i.e.,  parallel 
to  the  wavefront).  We  take  the  coordinates  of  these  receivers  to  be  (L,0) 
and(L,l).  By  the  same  procedure  as  was  used  to  derive  equations  (29)-(30) 
and  (74)- (75),  we  see  that  for  this  case 


f  L  L  L  A.l^'-C^^A^L'-  \')  NlrOJi-n.'dn'  j* /  (92) 

JT-.  --  f  I.L  L"  L„  L,  a»U'-  s')  m‘)w 


(93) 


where  £’  =  KJl. 

Following  the  same  procedure  as  in  equation  (29) -(41)  we  obtain 


S,^  =  t  f  (I.  +1.  ) 

BTfc; 

where  now 


i,=  CU' 

S.«  A,UC-iI\/wS')N(r')Jh  .h;  dr.' 

(96) 

vc'tr 

C  a.  e:  -  ^  ^  ♦  n  nu')  ^  a  t:  dc ; . 

(97) 

(94) 

(95) 


In  order  to  evaluate  lj_,  we  make  the  substitutions  (43)  and  (44)  and 
allow  the  integration  over  £  to  extend  from  -®  to  +03,  whence 


i.*  A.d^i^r) Nir'jd^iun 

r  IL'  L  Lo„  I')  N(r')  Ji| 


run1 


<98> 
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The  second  integral  in  (.98)  may  be  written  as 

j  -n  \  ,oo  i  ,  j_  i_  s  >^2*  -ri 

av-y,  a.?,  +  h  y  —  <&yy>  ) - IV  1“  o' 1  t  i f, )  +  e,  +  i  ( ^  "7 ' J  J  ri 


=  <&<■  ^  n^A  l1  L,  +•  * )]  L„  e  *"  J-r,  } 


(99) 


where  p  =  — ,o  -  4r  ,  ii1 

a  it,  and  q  =  — - — .  Since  the  real  part  of  p  is  positive,  we 
may  write  & 

-  ( y  V  )  -i—  -  “  '  r-  —  _^£—  '■  * 


s: 


^ , 


J—  f «  -  (  Tv  -n  -  -*pr  )  • 

hp)  e  J-  rP  ’ 


Jn 


HP 


\  ^  ^  \ 

J-c  Kr 


-  A_ 

‘t  £  ^ 


-  5-’ 


I  - _  A  ; 

H  0  J5,  1 


'  'J  .  ?, 


Substituting  (100)  into  (9-9) ,  we  obtain 

r  “■  -  pt  .  /  Vvri~.  -r_  \  . —  , _ i _ _ 

)  -  «,  ^  "  o*-"  V  li  +  vi1!  -  -i#  In  \ 


/Ji 


*^1 


‘  v  %  +■ 


(100) 

r-/iy  -j 

'J  r,U"-k  J  / 


r —  f  /  -  I  *•  h  V  j  t'  V  a'’*’  ) 

=  jn  JU)  Jt/V-Wat,  e  j) 


(101) 


and  our  expression  for  1^  becomes 


,  a,.  <  \_l1AIcV_1  -  ^ 

I,-;  1L  -Me*  -  I.  *nLu/.'M  Je  At} 


(102) 


It  is  not  apparent  from  the  form  of  equation  (102)  that  the  expression 
for  1^  reduces  to  the  value  found  for  the  single-receiver  case,  equation  (54), 
in  the  limit  as  the  separation  between  the  receivers,  l,  approaches  zero.  In 
order  to  make  the  relation  between  (102)  and  (54)  more  apparent,  we  substitute 
2£ 

C,  =  — 2  and  remove  the  complex  numbers  from  the  denominators  in  (102),  obtaining 
3 


T  ^  r  \  u?-/a'— > 7^1]  -  > 

V-L«  L\e  )„  stt;  nwj.1  Je  h} 


(103) 


We  mav  set 


(  +  »■  '  r  (  Off,  0  +  Jk  CVr  ft  ) 


(104) 
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r" 


where  r  *  +  1,  cos  ©  «  ^  ,  and  sin  ©  **  De  Moivre’s  theorem  then  implies  that 


^  ♦  1  -  ir  ( 


and  hence 


U>.  1  +  ).  jtk-V*'  1 


) 


< —  /  [  ■  i -  'i 

-  * r  U  ^T~  +  *  -j  x  / 

+c)l"x  +ilt  lIcTT  -OV . 


(105) 


Jr-*V.  1  =  t(ic\l  +Ol  +r(Kx  +  l+<'Jt  + 


-  x  -1<*£  +/31  -t-  j[  -J  <*“■  +2«/3 

where  a  =  (/?  +  1  +  and  6  =  (Z^2  +  1  -  5)^".  This  last  expression 


(106) 


reduces  to 


-5T7*  t’1  11  =  i  [1  'i)1'  +  i  (iIc\T  +  } 

so  equation  (103)  becomes 


(107) 


L'a'X  f  ‘ 

I,-'  "IT  ). 


i  C  •  + 1 


8  5 ^  —  i.  ►  n/l/  \ 

e‘  e  ’ 


(108) 


which  reduces  to  (54)  when  !'  =0. 


a 1 2g2 


If  the  restriction  imposed  by  (17)  is  valid,  the  term  e  4  will  cause 
the  integrand  to  vanish  for  £  >1,  so  throughout  the  effective  range  of  inte¬ 
gration  we  may  neglect  £2  im  comparison  with  unity.  For  this  case  equation 
(10 8)  reduces  to 


I,'-  LV"  c 

FT  We  c"  “ 


°'v  J 

I'1  8 


-  - - *■  r  t'1  v 

e  11 


r  SL 

r  .  %.  o.v 

-  La k  e  e 


(109) 

which  is  the  same  value  as  that  which  Chernov  derived  for  the  three-dimensional 
case  (p.  102,  (213)).  For  small-scale  inhomogeneities  (i.e.,  those  for  which 
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equation  (17)  is  not  valid)  we  shall  have  to  evaluate  (108)  numerically. 


We  may  evaluate  I2  by  substituting  (43)  and  (44)  into  (97)  and  setting 
the  limits  on  the  integral  over  £  equal  to  and  +*,  thereby  obtaining 


Ix  =  C  L-  A,  +  )  N<r ')  Jn  H  lx. 

Inspection  of  equations  (98)-(102)  enables  us  to  write 

l  \  -12 


(110) 


.  ■,  i  (iV*  >  _  \  1  >'  ! q  ~ _ _  i-l-  . 

1  "T  ( u'  I  r  -1 

=  1.  j  a*  { 


flr. 


i!  rrr 


an) 


where  we  have  used  the  definition  of  the  wave  parameter  D,  equation  (60). 
If  we  substitute  t  =  -  ii1 2/a* 2  into  equation  (111),  we  find 

4  -  i 

i  •.  =  ~ c  s’’  r!jt! 


where  T(a,e)  is  the  incomplete  gamma  function,  defined  by 


(112) 


e-1 


i 


t*  -  I 


. 


(113) 


We  now  show  that  equation  (112)  reduces  to  the  value  which  was  found  for  the 
single-receiver  case,  equation  (62),  as  the  separation  between  receivers,  i, 
approaches  zero.  Expanding  the  incomplete  gamma  function  in  a  power  series, 
we  find 


-  u.  (n-i)-  i 


(-O  l 


ni  (-i  +  r) 


} 


-ll  * 

> 

hence 


(114) 
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(115) 


fn  *?  ,  ,  -*il  I  *  /g*  1  ,  ,  r  L  \ 

1  1  (  D-i  '  ~  | 

|-<C 

^Yduni') 

In  *?•’  v  ,  i  k 

'  Up  tO  v'l] 

=  r^r  *io  i^)i(ui+»rv-i]^- 

which  agrees  with  (62).  Having  demonstrated  the  behavior  of  for  small  l, 
we  may  now  delete  the  first  term  inside  the  braces  in  equation  (112),  since 
it  has  no  imaginary  part  and  contributes  nothing. 

Since  the  incomplete  gamma  function  is  not  tabulated  for  complex 
arguments,  it  will  be  necessary  to  evaluate  I2  by  means  of  the  power  series 
expansion  which  is  shown  in  (114).  Because  this  series  is  only  slowly  con¬ 
vergent  when  the  wave  parameter  D  is  small,  we  will  find  it  helpful  to  treat 
this  particular  case  by  expanding  (112)  in  a  Taylor  series: 


(116) 


where 


= 


_  il 

D-i.  '  ex' 


.  n  il  -  nv  il  ■  rl  il 


(117) 


We  see  from  the  definition  (113)  that 


r,  /  .  “ 'h  oc  -  l 

V  («,*)  =  -  e  ? 


(118) 


r"(«,*)=  e'V'1  +(i-Oe'ir 


(119) 


Retaining  terms  of  up  to  third  order  in  the  small  quantity  D,  we  find 


lx*  o  5  ”  [D  1  +  I«  ’  si)].  !>< 


(120) 
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As  was  done  for  the  case  of  longitudinal  receiver  separation  in  equations 
(90)  and  (91),  we  write  the  expressions  for  the  correlation  coefficients  for 
phase  fluctuations  and  for  amplitude  fluctuations: 


s, trans 


*b, trans 


I,(o)  +  Ijo) 

1,(»  <  - 


(121) 


(122) 


We  note  that  we  may  obtain  convenient  expressions  for  the  correlation 
coefficients  in  the  particular  case  D  <<  1  by  expanding  (62)  in  a  power  series: 


\  _  kl  ii-  '] 

— —  -1  +  -  c  +  *  V  ...  -  1  ] 


--  ^  Kalall4  -A'  T  4  T  +  ...)-U'  TT  +•••’)+  ..  j 


ul  _X£—  1 

-  4T  yOl.aU  -  <B  +  >14  +...  J  . 


Substituting  (109),  (120),  (52),  and  (123)  into  (122),  we  find 
c  e~  e  6  ^  \_0  T  4  P*  (-  L?  -V  Ti:  ~  ft)  j 


~  /  --.r  4$-  ' 

■  (  -  'ar  4  1- 


D«  1  ,  Ka»l,  Ifca. 


(123) 


For  the  corresponding  three-dimension  case  Chernov  finds  the  value  (p.  102, 
(227))  Rb  =  (n‘*/2a‘*  -  2£2/a2  +  l)exp(-£2/a2) .  Similarly,  equation  (121)  in 
this  particular  case  reduces  to 


~  e 


(125) 


which  is  the  same  value  found  by  Chernov  (p.  107,  (227a)). 


Summary  of  Theoretical  Results 


The  expression  for  B2,  S2,  p,  Rb,long(AL>»  Rb,  trans (Z)  >  and  Rs,trans(*} 
are  summarized  arid  compared  with  Chernov's  values  for  tne  three-dimensional 
case  in  Table  I  of  Rivers  and  von  Seggern  Q97JJ)  .  It  may  be  seen  that  the  ratio 
of  the  amplitude  and  phase  fluctuations,  BVS2,  and  the  correlation  between  them, 
p,  can  be  expressed  parametrically  as 


B2 


g2  =  fi(K.L.a) 

(126) 

p  =  g1(K,L,a). 

(127) 

If  we  make  the  approximation  Ka  >>  1,  these  equations  take  on  the  even 
simpler  form 

rr  =  f,(D)  (128) 

S2  2 


P  =  g2 (D) •  (129) 

Figures  2,  3,  and  4  show  B2,  S7,  and  p  as  functions  of  D.  We  may 
eliminate  D  from  (128)  and  (129)  and  plot  B^/S^"  versus  p,  in  Figure  5. 

The  values  of  D  corresponding  to  selected  points  along  this  curve  are 
shown  on  it.  Analogous  plots  may  be  obtained  from  the  more  exact  rela¬ 
tions  (126)  and  (127)  by  taking  K  and  L  to  be  given  and  eliminating  a.  This 
is  done,  for  three  different  source-to-receiver  distances  L,  in  Figu£es_6, 

7,  and  8.  Comparison  of  these  three  plots  shows  that  the  values  of  B2/S2  and 
p  are  insensitive  to  the  value  of  L,  so  long  as  L  >>  a.  It  was  this  insensi¬ 
tivity  to  the  length  of  the  path  through  the  scattering  medium  (after  the 
incident  wave  has  been  scattered  several  times)  which  enabled  us  to  assume 
that  all  paths  were  of  the  same  length  L  when  computing  the  averages  (29) 
and  (30). 


By  once  again  assuming  K  and  L  to  be  given  constants,  we  may  calculate 
the  longitudinal  amplitude  and  phase  correlation  coefficient,  R^(AL)  and  RS(AL), 
as  functions  of  a.  This  is  done  in  Tables  I  and  II  for  separation  between 
receivers  of  up  to  100  km,  approximately  the  distance  across  NORSAR.  It  is 
seen  that  both  R^(AL)  and  Rg(AL)  fall  off  slowly  with  increasing  separation 
between  receivers.  (The  apparent  increase  in  Rfc  with  increasing  AL  for  large 
values  of  a  is  an  artifact  of  the  computation  process;  in  fact,  R^(AL)  z  1.0 
for  AL  <<  a.)  The  tabulated  values  were  computed  for  source-to-receiver 
distance  L  =  3000  km;  they  were  found  to  be  scarcely  different  for  the  case 
L  =  9000  km. 

The  transverse  amplitude  and  phase  correlation  coefficients,  Rb(£)  and 
Rs(£),  were  similarly  calculated  for  fixed  K  and  L  and  for  various  values  of  a. 
As  shown  in  Tables  III  and  IV,  Rs(Jl)  fall  off  in  an  approximately  Gaussian 
curve  with  increasing  separation  between  receivers,  but  ( JL )  decreases  to 
negative  values  for  l  z  a.  Inspection  of  Table  V  shows  that  the  values  of 
Rb ( ^ )  are  relatively  insensitive  to  the  assumed  value  of  L;  Rs(£)  is  even  less 
sensitive  to  L.  The  behavior  of  R^(f)  and  of  Rs(£)  are  shown  graphically,  for 
both  large  and  moderate  values  of  the  wave  parameter,  in  Figures  9-12.  These 
graphs  demonstrate  the  same  general  behavior  with  increasing  9,  as  do  those 
given  by  Chernov  for  the  three-dimensional  case.  (Chernov,  Figs.  5  and  6) 

Figure  13  shows  R[,(S,)  for  five  values  of  D  and  extends  the  separation  SL  to 
distances  beyond  the  100  km  considered  in  Table  III. 
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Figure  2.  Theoretical  logarithmic  amplitude  fluctuation  as  a  function  of 
wave  parameter  D. 


3-Dimensional 


Theoretical  correlation  coefficient  of  logarithmic  am] 
fluctuations  and  phase  fluctuations  (at  a  single  rece 
function  of  wave  parameter  D. 


AMPLITUDE  FLUCTUATION  TO  PHASE  FLUCTUATION 


FLUCTUATION  RATIO  B**2/S**2 


Figure  6.  Theoretical  values  of  — j  and  p  for  various  characteristic  dimensions  A  of  the  random 

......  S 


scatterers. 
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FLUCTUATION  RATIO  B**2/S**2 
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TRANSVERSE  SEPARATION  BETWEEN  RECEIVERS  /  IN  UNITS 
OF  CHARACTERISTIC  DIMENSION  A  OF  SCATTERERS 


Figure  9.  Theoretical  transverse  amplitude  correlation  coefficient  for  case  D<<1. 


TRANSVERSE  SEPARATION  BETWEEN  RECEIVERS  /  IN  UNITS  OF 
CHARACTERISTIC  DIMENSION  A  OF  SCATTERERS 

Figure  10.  Theoretical  transverse  amplitude  correlation  coefficient  for  case  D  »  10. 
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TRANSVERSE  RS(D<<1) 


TRANSVERSE  SEPARATION  BETWEEN  RECEIVERS  /  IN  UNITS 
OF  CHARACTERISTIC  DIMENSION  A  OF  SCATTERERS 

Figure  11.  Theoretical  transverse  phase  correlation  coe:  f  ii  icnt  K  ( 1 )  for  cast-  D<<1. 

s 


AI-38 


TRANSVERSE  R,  (0=10) 


TRANSVERSE  SEPARATION  /  BETWEEN  RECEIVERS  IN  UNITS  OF 
CHARACTERISTIC  DIMENSION  A  OF  SCATTERERS 

Figure  12.  Theoretical  transverse  phase  correlation  coefficient  Rg(l)  for  case  D  =  10. 
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APPENDIX  II 

EVALUATION  OF  CERTAIN  INTEGRALS  USED  IN  APPENDIX  I 


,  t 


— 1 — ^ 1 1  pin  ml  *  c-**-  -r  ~ 


APPENDIX  II 


Evaluation  of  certain  integrals  used  in  Appendix  I 
In  this  appendix  we  evaluate  three  integrals  involving  the  functions 


A(  (<*, - 


(A.  1) 


and 


t  j  f--x  2—  I 

+  1/ 

which  appear  in  the  two-dimensional  modification  of  Chernov's  theory, 
three  integrals  are 


(A. 2) 
The 


Aama,  =  )_«,  A,  (*, ,'j,+  p  A 1  Ui^  +  y)  Jy 


and 


1 


r  00 


It  follows  from  the  definition  (A.l)  that 


-W, 


-1 


OO  A  8,  -‘®I 

e  -  e 

>6, 

e  -  e 

kk 

kk 

« IV  8,  > 

e 

-a„a('  v'r-k'*'  j-~  iv 


T  f"  +  „  r0'  ....  , 

s  jy-^u  e 


-  e  +  e 


J  <*y 


where  we  have  made  the  substitutions 


ft  -  iiilT  .  JL 

“1  “  IK,  +  4 


AII-1 


(A.  3) 


(A. 4) 


(A.  5) 


(A. 6) 


(A.  7) 


VvV., 


and 


I Ji+I) 1  ,  JL 

M  +  *• 


(A. 8) 


In  order  to  evaluate  the  two  integrals  appearing  in  (A. 6),  we  make  use  of 
the  identity  .i 


V.  (»  -cbji  +  C  ) 


=Me'a>r/^4 

=  K.(e‘,c'fe)J¥  £*'} 


3L 
1  a 


I  _kL  5L  \  -  . 

-  (.£  -  4a  +  H  /  .  o  >0  j 


expressions  appearing  in  the  exponents  in  (A. 6)  are 


0.*0  -  t 


[  i  ( J-  .  ±  )  l  IK  Hi  \  1L  lx  1  r 

M  0,1  +/'-ot>  +  01 1  +  J  +  x 


(A.  9) 


(A. 10) 


jt  .  /  _l.  _  _i_  \  (Hr  _  Hi.\  lx  _  lx.  1 

0,-0t'*Lyl<*,  <*-J+yw,  oij-t  k,  '‘i  J , 

so  we  may  evaluate  the  first  integral  in  (A. 6)  by  substituting 


(A. 11) 


a  -  'a<n,o <t 

(A. 12) 

uMi+ol,/ 

i 

«.y>  +«»/ 

,*  3_ 

(A. 13) 

+  5. 

(A. 14) 

into  (A. 9).  It 

follows  from  the  definitions  (A.12)-(A.14)  that 

b1  .  2_ 

(y.-y)1 

C  "  ho  +  1 

1  (<*,+«  jJ  4  ; 

(A. 15) 

AII-2 


so  we  have 


;ifr. 


0*1. 


a(ot,+«,) 


+ 


-  _  { 

J  «.♦*»  ""Wfci, 


(A. 16) 


In  order  to  evaluate  the  second  integral  in  (A. 6),  we  note  that  the 
restriction  a  >  0  was  used  in  deriving  (A. 9),  so  we  may  substitute  from 
(A. 11)  into  that  identity  only  if  ot2  >  aj.  We  may  treat  the  case  04  >  a2, 
however,  simply  by  noting  that 


<*  <» 


(A. 17) 


so  we  may  treat  both  the  cases  a2  >  aj  and  04  >  a2  by  substituting 


a  - 


W,  -*-l\ 
lo(  *01-, 


(A. 18) 


into  (A. 9).  We  find  that 

C  ”  H4  +  1  ■  }  \o<1-  «1  J  +  *1 

and  thus 

- 

f.  J_*  e  A 


rSP  id-hllt 


j  ZlLSlil 


+  t) 


(A. 19) 

(A. 20) 

(A. 21) 

(A. 22) 
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Substituting  (A. 16)  and  (A. 22)  into  (A. 6),  we  derive  as  our  first  result 


_  _ . —  .  I  j •  I- 

Ia,/,'  p.  [7r<«,Wi  V  +  lh-M*,-*,!  V  +  *•  / 


=  +  A  i  )1 . 

We  may  use  (A. 16)  and  (A. 22)  to  evaluate  (A. 4): 

.  ,  ^  te,  -i»t 

-r  _  . _ * -  f  e  *  e _  e  *  e _ 

J',.  i^'oT  X  3l  / 

_  _ i _  p  T  ^'O  ,l-e,-ol 

le  +  e  +e  +  e  J, 


(A. 23) 


'  V  ( “ 
h-.  PP  1 **  ^ 


Kl.e  +  Fu,  S _ e*1  ‘ 


;  L  A,Ul  +  at>y1-'/i^  +  A,  >),')■.)]  , 


For  CA.5)  we  find 


T  _  _ L —  p  5 - - 

A i 2 T  Jrf,ol t  ) - oo  1 


‘9,  -i*.  p',6>- 

e  1  -  e  e  -  e  ... 


Jl  d) 


- p  \  iUt-e,)  ;l-vOl 

va  ia.tr  >_<>.  \.e  +  e  -  e  -  e  J 

- — i=W  V“  P‘l(9'+®^  .  c“  1 

-  4*  !<*,*,.  e  Ay  +  ^  e  Jiy  J, 


By  inspection  of  (A.  9.)  we  may  write 


*>  ii«yv*v>>|*c)  nr  .  ,  v_  j5.\ 

4kv,  )_ „  e  J[y  -  H  a  /OV*-  VC'Ma  +1/  i^a^>Ot 


(A. 24) 


(A. 25) 


(A. 26) 


Our  derivation  proceeds  as  before,  but  now  our  treatment  of  the  case  aj  >  a2 
replaces  (A. 17)  with  the  identity 


M  c  A  v  =  ~  „  6  A  > 


(A. 27) 


and  our  final  result  is 


AII-4 


(A. 28) 


a»"  +  -  AviU.-cxJj  7,-y,)] 

where  the  plus  sign  Is  to  be  used  If  a?  >  otj,  and  where  the  minus  sign  is  to  be 
used  if  ai  >03. 
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